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ABSTRACT 


N-iEarth  cover  provides  protection  to  underground  structures  against 
the  effects  of  air  shock  loading.  Part  of  this  protection  may  result 
from  attenuation  of  stress  with  thickness  of  the  cover.  Measurements 
of  vertical  earth  stress  at  three  depths  and  at  five  ground  ranges 
were  made  during  Shots  9  and  10  of  Operation  UPSHOT -KNOTHOLE  to  detect 
and  evaluate  stress  attenuation  with  depth.  ."’Data  fit  equally  well  the 
empirical  equations  ,  \  \ 


?  =  P,  exu-  (d  -  d. )  0.07 
1  *  \  •  X' 


P  =  P, 


/  d_\  -°-3^ 


in  which  P  add  P]_  are  the^ stresses  in  psi  at  depC^s  d  and  dj_  in\ft,  d^_ 
Being  the  shallower-  Precision  in  each  -case  is  better  than  +  25\per 
cent.  \  \  \  \  \ 


Fart  II  -  Project"  1.4b 


A  practical  system  for  measuring  free-field  earth  stresses  and 
strains  resulting  from  transient  loads  has  been  tested  with  sufficient 
thoroughness  to  establish  its  feasibility.  Arrays  of  directionally 
sensitive  earth  stress  and  strain  gages  and  accelerometers  were  in¬ 
stalled  5  ft  deep  to  record  these  paramet  rrs  during  Shots  1,  9>  and  10 
of  UPSHOT -KNOTHOLE.  Duplicate  instrumentation  showed  tnat  stress  meas¬ 
urements  were  reproducible  with  average  deviations  of  l6  per_cent  and.— -s 
strain  measurements  with  average  deviations  of  35  per  cent This  teht 
disregards  the  presently  unknown  factors  related  to\ pertnr nttions  of  \ 
the  stress  field  bj\the  gages.  Stress-strain  graphsX demonatkeatV  hys-  \ 
teresis.  'plastic  deformation  of  the  soil\resulted  ir\  rates \of^  energy 
dissipation -as  high  at  300  pft-lb/ft3.  Da\a  from  ShotvlO  de^dXed  Yhe 
stress  tensor  in  terms\of  magnitudes  and  directions  of  \he /thr^  prin-, 
cipal  stresses  as  a  function  of  time.  \  \/  \ 
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FOREWORD 


Measurement  of  free-field  phenomena  in  the  earth  resulting  from 
airbursts  of  nuclear-  weapons  was  divided  under  Project  1.4  of  UPSHOT- 
KNOTHOLE  into  two  subprojects:  Project  1.4a  was  concerned  with  the 
detection  and  evaluation  of  attenuation  of  vertical  earth  stress  with 
depth  and  was  intended  to  give  information  pertinent  to  Project  3*8 
structural  response  tests;  Project  1.4b  involved  testing  instruments 
and  analytical  procedures  for  determination  of  stress  and  strain  tensors 
produced  within  the  earth  as  a  result  of  a  high-yield  explosion. 

Part  I  of  this  report  deals  entirely  with  Project  1.4a  measurements 
of  vertical  stress  as  a  function  of  depth,  time,  ground  range,  and  inci¬ 
dent  air  overpressure.  Part  II  is  concerned  with  Project  1.4b  systems 
for  measurement  and  analysis  of  data  for  determining  stress  and  strain 
tensors  in  the  earth  as  a  function  of  time. 

This  is  one  of  the  reports  presenting  the  results  of  the  70  pro¬ 
jects  participating  in  the  Military  Effects  Tests  Program  of  UPSHOT- 
KNOTHOLE,  which  Included  11  test  detonations.  Reference  to  other  per¬ 
tinent  test  information  Is  made  in  WT-782,  Summary  Report  of  the 
Technical  Director,  Military  Effects  Program,  which  Includes  the  fol- 

T - -t - - 

a.  An  over-all  description  of  each  detonat_on  including  yield, 
height  of  burst,  ground  zero  location,  time  of  detonation, 
ambient  atmospheric  conditions  at  detonation,  etc.,  for  the 
11  shots . 

b.  Compilation  and  correlation  of  all  project  results  on  the 
basic  measurements  of  blast  and  shock,  thermal  radiation, 
and  nuclear  radiation. 

c.  Compilation  and  correlation  of  the  various  project  results 
on  weapons  effects. 

d.  A  summary  of  each  project,  including  objectives  and  results. 

e.  A  complete  listing  of  all  reports  covering  the  Military 
Effects  Tests  Program. 
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PART  I 

VERTICAL  STRESS  AS  A  FUNCTION  OF  DEPTH  AND  RANGE 


1.1  HISTORICAL  BACKGROUND 


Experience  has  shown  that  buried  structures  are  damaged  less  by 
air  blast  than  are  similar  uncovered  structures.  Navy  structure 
3.2.6  1/survived  Easy  shot  of  Operation  GREENHOUSE  essentially  un¬ 
damaged  when  covered  with  loose  sand  to  a  minimum  depth  of  2  feet,  but 
was  severely  damaged  by  air  shock  of  approximately  the  same  intensity 
on  Item  shot  after  the  earth  cover  had  been  removed. 

The  mechanics  of  the  protective  effect  of  earth  cover  is  not  com¬ 
pletely  understood.  Protection  afforded  by  minor  thicknesses  of  earth 
is  almost  entirely  a  result  of  altered  structural  response  since  the 
earth  cover  increases  the  effective  mass  of  the  structure,  decreasing 
its  response  to  transient  loads .  Stability  of  earth  cover  over  struc¬ 
tures  which  extend  above  the  ground  surface  requires  side  slopes  of 
30-45°.  Such  slopes  have  been  found  to  reduce  the  reflection  inten¬ 
sity  of  air  shock  waves  to  about  50-75  per  cent  of  that  caused  by 
vertical  walls. 

Much  of  the  protection  afforded  by  earth  cover  over  deep  under¬ 
ground  structures  is  derived  from  radically  different  phenomena.  It 
may  result  In  part  from  attenuation  of  stresses  within  the  soil  column 
and  in  part  from  development  of  arching  above  the  structure .  The 
latter  condition  results  when,  as  a  consequence  of  loading,  a  portion 
of  the  structure  such  as  the  roof  slabs  yield  sufficiently  to  cause 
shear  stresses  to  develop  within  a  zone  of  soil  overlying  the  struc¬ 
ture  and  extending  laterally  beyond  it.  The  shear  stresses  thus  de¬ 
veloped  divert  part  of  the  vertical  load  from  the  roof  to  the  sur¬ 
rounding  soil.  This  process  acts  as  though  arches  had  developed  in 
the  soil  above  the  yielding  member  of  the  structure  to  transmit  the 
load  to  adjacent  soil  or  parts  of  the  structure.  Occurrence  of  arch¬ 
ing  depends  primarily  on  a  sufficient  thickness  of  soil  of  the  proper 
type  to  support  relatively  high  shear  stresses.  Movement  of  the  struc¬ 
ture,  particvjlarly  of  a  vibratory  nature,  or  plastic  flow  of  the  soil 
will  decreaae  the  stability  of  arching  and  reduce  its  effectiveness. 

The  thickness  of  cover  which  will  distinguish  between  these 
modes  of  protection,  mass  effect  or  arching,  is  vegue  and  debatable. 
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It  is  safe  to  anticipate  no  effective  arching  for  earth  cover  shallower 
than  the  minimum  lateral  dimension  of  the  covered  structures.  However, 
there  do  not  appear  to  he  any  data  which  indicate  the  relation  between 
the  span  of  the  structure  and  thickness  of  earth  cover  for  which  soil 
arches  will  develop  under  dynamic  loading.  Moreover,  it  seems  probable 
that  if  earth  stresses  are  attenuated  with  depth,  effective  attenuation 
may  be  anticipated  for  any  thickness  of  cover  sufficient  to  offer  pro¬ 
tection  through  arching  over  a  structure. 

Attenuation  of  stress  and  arching  depend  upon  mechanical  proper¬ 
ties  of  the  soil.  The  nature  of  soils  suggests  that  if  there  is  rela¬ 
tive  movement  between  adjacent  soil  elements,  dissipation  of  energy 
and  consequently  attenuation  of  stresses  may  be  anticipated.  Theoret¬ 
ical  studies  of  earth  stress  and  strain  phenomena  resulting  from  steep- 
fronted  transient  loads  of  several  hundred  milliseconds  duration  are 
not  yet  well  advanced.  Some  doubt  has  been  expressed  that  any  de¬ 
crease  in  stress  with  depths  of  the  order  of  5-15  ft  should  be  antici¬ 
pated  . 

Measurements  of  vertical  stress  at  depths  a3  great  as  20  ft  were 
made  during  Operations  BUSTER -JANGLE  and  TUMBLER-SNAPPER  to  determine 
experimentally  whether  attenuation  of  air-shock-induced  stress  could 
be  observed.  Results  from  BUSTER -JANGLES/  were  so  erratic  that  no 
conclusive  evidence  of  attenuation  could  be  deduced.  Inadequate  con¬ 
trol  of  backfill  compaction  appeared  to  be  largely  responsible  for  the 
erratic  results. 

Gage  arrays  of  two  different  types  were  used  to  measure  vertical, 
earth  stress  at  several  depths  during  TUMBLER -SNAPPER .1/  Performance 
of  gages  placed  in  individual  holes  was  unsatisfactory  because  of  un¬ 
solved  placement  problems.  Other  gages  placed  in  a  mound  of  compacted 
earth  apparently  performed  satisfactorily,  but  unanticipated  oscilla¬ 
tion  of  the  mound  obscured  conclusive  evidence  that  attenuation  of 
vertical  stress  occurred. 

Free-fleld  measurements  of  another  parameter,  vertical  component 
of  acceleration,—'  have  indicated  a  decrease  in  peak  acceleration  with 
depth  although  the  data  are  too  limited  in  number  and  scope  to  do  more 
than  indicate  a  trend. 

The  status  of  the  problem  of  measuring  attenuation  of  air-shock- 
induced  earth  stress  with  depth  was  practically  unaltered  between 
BUSTER- JANGLE  and  planning  for  Operation  UPSHOT -KNOTHOLE.  However, 
both  previous  operations  had  indicated  the  critical  importance  of 
proper  compaction  of  backfill  around  and  above  the  gages. 

1.2  PLAN  AND  PURPOSE  OF  PROJECT  1.4a 


Measurement  of  the  vertical  component  of  earth  stress  at  several 
depths  and  ground  ranges  was  included  in  plans  for  UPSHOT -KNOTHOLE  as 
Project  1.4a.  These  measurements  were  designed  primarily  to  determine 
whether,  with  adequate  control  of  backfill,  attenuation  of  vertical 
stress  could  be  observed  at  depths  of  15  ft  or  less.  Identical  gage 
arrays  at  five  ground  ranges  were  planned  to  improve  chances  for  ob¬ 
serving  attenuation  and  indicate  possible  dependence  of  such  attenua¬ 
tion  on  the  magnitude  of  the  locally  incident  air  shock  wave.  Location 
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of  the  arrays  was  planned  to  Include  anticipated  incident  peak  over¬ 
pressures  ranging  from  80  to  10  psi . 

Data  from  these  measurements  were  expected  to  be  comparable  with 
concurrent  free-field  measurement  of  vertical  stress  at  several  depths 
in  the  earth  fill  adjacent  to  the  buried  structures  of  Project  3*8*  5 / 
It  was  anticipated  that  such  comparison  might  provide  a  correlation  be¬ 
tween  free-field  earth  stresses  and  structural  damage. 

1.3  INSTRUMENTATION 


Carlson-Wiancko  earth  stress  gages  of  improved  design,  described 
in  Part  II  of  this  report.,  were  used  for  all  but  three  measurements  on 
Project  1.4a;  older-type  gages  were  used  in  these  three  instances  to 
determine  whether  there  was  any  radical  difference  in  performance  of 
the  two  types.  Response  of  each  gage  was  transmitted  over  buried  cable 
to  a  recording  shelter,  where  it  was  recorded  by  a  Consolidated  Engi¬ 
neering  Corporation  carrier-amplifier  photographic  recording  system 
backed  up  by  Ampex  magnetic  tape  recorders .§/ 

Gages  were  placed  at  depths  of  1,  5j  and  15  ft  below  ground  level 
at  distances  of  750,  1000,  1250,  1420,  and  2000  ft  from  nominal  ground 
zero  for  Shots  9  and  10  (Fig.  l.l). 


* .  ..  rl 

- I — 1 -  - . 

F-203  F-282 


F-21B 


r*  -  2500 


.  -  '  j  1  ft  deep 

•  - 5  ft  deop 

- 10  ft  docp 


Nominal 
ground  zero 


Fig,  1.1  Instrumentation  for  Earth  Stress  Attenuation  -  Shots  9  and  10 
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These  distances  included  ranges  of  predicted  peak  air  overpressures 
between  25  and  12  psi  from  Shot  9  and  between  80  and  10  pal  frost  Shot 
10.  Incident  air  overpressures  at  ground  level  were  obtained  from 
measurements  by  Project  1.1b  and  from  a  pair  of  ground-baffle  gages  at 
Station  F-28l  for  Project  1.4b. 

Holes  for  the  earth  stress  gages  were  large  enough  to  permit 
direct  hand  tamping  of  the  backfill  Immediately  surrounding  each  in¬ 
strument.  Borings  for  the  15-ft  gages  were  of  sufficient  diameter  to 
permit  use  of  30- in.  diameter  corrugated  iron  casing  to  make  the  holes 
safe  for  occupancy  by  personnel. 

All  gages  were  carefully  bedded  and  leveled  in  the  bottoms  of 
the  holes,  which  had  been  previously  tamped  and  trimmed  to  within  about 
0.1  ft  of  the  prescribed  elevation.  The  final  elevation  of  each  gage 
was  determined  with  a  precision  of  0.01  ft,  and  the  spread  of  eleva¬ 
tions  for  all  gages  at  a  particular  nominal  depth  was  held  to  0.24  ft. 
Screened  moist  soil  was  hand  tamped  around  the  gages,  special  precau¬ 
tions  being  taken  to  ensure  uniformity  of  the  fill  and  to  bond  the  fill 
to  the  relatively  undisturbed  soil  surrounding  the  hole  by  wetting  the 
sides  of  the  hole  immediately  prior  to  placing  the  fill. 

In  the  deep  holes  the  initial  13-ft  length  of  casing  was  with¬ 
drawn  and  replaced  by  shorter  lengths  aa  backfilling  progressed.  A 
large  pneumatic  tamper  was  used  after  about  12  to  l8  in.  of  backfill 
had  been  hand  tamped  over  the  gages  to  provide  reasonable  protection. 
Backfill  was  placed  and  tamped  in  lifts  of  about  6  In.  Uniformity  of 
the  compacted  backfill  to  a  height  of  somewhat  more  than  one  hole  di¬ 
ameter  was  regarded  as  extremely  important,  but  as  a  rule  reasonable 
precautions  were  taken  to  maintain  uniformity  of  compaction  throughout 
the  entire  column.  Placement  of  the  gages  at  the  1-  and  S-fb  depths 
differed  from  that  of  the  deeper  ones  only  in  that  the  sides  of  the 
holes  were  sloped  to  minimize  arching  over  the  gagas.  Backfill  was 
built  up  to  a  height  of  about  6  in.  above  nominal  ground  level  and  cut 
back  to  grade  during  final  preparation  of  the  blast  line  area. 

All  backfill  was  prepared  by  blade-mixing  local  surface  soil  with 
water  to  make  it  somewhat  wetter  than  optimum  for  compaction  because 
the  dry  windy  climate  and  blowing  dust  tended  to  dry  the  material  very 
rapidly.  As  placed  and  tamped  it  was  probably  on  the  dry  side  of  opti¬ 
mum.  Compaction  of  drier  than  optimum  backfill  was  desirable  to  pre¬ 
vent  separation  of  backfill  from  the  undisturbed  soil  as  well  as  to 
prevent  shrinkage  away  from  the  gages. 

A  few  random  soil  density  measurements  Indicated  that  the  com¬ 
pacted  fill  was  about  10  lb/cu  ft  more  dense  than  the  undisturbed  soil. 
No  measurements  were  made  to  determine  relative  compressibilities  or 
moduli  of  undisturbed  and  compacted  soil.  Presence  of  soluble  salts 
in  the  playa  soil  of  Frenchman  Flat  made  the  low-density  undisturbed 
material  relatively  hard  and  rigid  as  a  result  of  cementation.  The 
same  factor  probably  caused  the  denser  compacted  fill  to  become  some¬ 
what  more  rigid  as  recementation  took  place,  and  it  is  possible  that 
the  compressibilities  may  have  been  more  closely  matched  than  the 
measured  densities  suggest.  Tests  performed  for  Project  3 .85/  show 
that  the  shearing  strength  of  compacted  soil  from  Frenchman  Flat  in¬ 
creased  by  a  factor  of  5  to  8  as  it  dried  from  optimum  water  content 
to  about  l/5  of  optimum. 
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1.4  PRESENTATION  OF  RESULTS 

Valid  records  were  obtained  for  all  Project,  1. 4a  gages  from 
Shots  9  and  10.  Copies  of  the  stress -time  curves  from  these  recoi'ds 
are  presented  in  Figs.  A.l  through  A. 5  of  Appendix  A.  Ground  and 
slant  ranges  for  each  station  shown  on  the  plots  are  approximate  true 
ranges  to  the  center  of  each  array  and  differ  from  true  range  to  a 
specific  gage  by  a  maximum  of  11  ft  for  ground  ranges  and  2  ft  for 
slant  ranges.  True  ranges  diverge  widely  from  the  nominal  ranges  used 
in  planning  Project  1.4a  because  of  the  large  error  in  ground  zero  for 
Shot  9  and  the  direction  of  the  error  for  Shot  10  (Fig.  1.2). 


Ac ( uni  ground  zero.  ,%ol  # 


Fig.  1.2  Relative  Positions,  Nominal  and  True  Ground  Zeros  -  Shots 
9  and  10 


Pertinent  data  from  these  gages  are  presented  in  Tables  1 .1  (Shot,  9) 
and  1.2  (Shot  10 ).  Incident  air  overpressures  measured  by  Stanford 
Research  Institute  for  Project  1.1b 1/  are  included,  as  are  those  meas¬ 
ured  at  Station  F-28l  as  a  part  of  Project  1.4b. 

Peak  earth  stresses  and  incident  air  overpressures  are  plotted 
as  functions  of  ground  range  in  Fig.  1.3  f°r  Shot  9  and  Fig.  1.4  for 
Shot  10.  These  curves  emphasize  coherence  of  the  air  overpressures 
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Fig,  13  Vertical  Earth  Stress  -  Fig,  1.4  Vertical  Earth  Stress  - 
Shot  9  Shot  10 


and  of  the  earth  stresses  observed  at  the  15-ft  gages.  This  coherence 
contrasts  notably  with  the  erratic  character  of  some  of  the  data  from 
the  1-  and  5-ft  gages.  Curves  were  drawn  through  all  data  points  ex¬ 
cept  where  deviation  from  a  smooth  curve  is  small.  Light  dashed  curves 
have  been  drawn  for  the  1-ft  data  from  Shot  9  and  the  1-  and  5~ft  data 
from  Shot  10  to  indicate  a  possible  smooth  curve  which  conforms  in 
general  to  the  shapes  of  the  incident  air  overpressure  and  the  15-ft 
gage  curves . 

1.5  ANALYSIS  AND  DISCUSSION  OF  RESULTS 

Data  from  the  earth  stress  gages  at  different  depths  will  be  use¬ 
ful  for  determining  variations  of  air-shock- induced  stresses  with  in¬ 
creasing  depth  if  it  can  be  established  that  the  major  source  of  these 
stresses  is  the  air  overpressure  incident  at  the  adjacent  ground  Bur- 
face.  The  criteria  of  adequacy  of  data  are  comparisons  of  arrival 
times,  durations,  and  stress  magnitudes  as  well  as  similarities  of  air 
overpressure-time  and  earth  stress-time  wave  forms. 
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Direct  comparison  of  air  overpressure  and.  earth  stress  data  is 
valid  only  for  comparable  ground  ranges .  Consequently  comparison  of 
air  -and  ground  data  for  Shot  9  (Table  l.l)  is  limited  to  those  from 
Jtations  F-28l  (air  only),  i>’-2 82,  and  F-200;  comparison  or  data  from 
Station  F-2l6  is  of  doubtful  value  because  of  the  large  difference  in 
ground  ranges  to  the  air  pressure  and  earth  stress  gages.  Similarly 
air  and  ground  data  for  Shot  10  (Table  1.2)  are  comparable  at  Stations 
F-2l6  and  F-200  but  are  of  doubtful  comparability  at  Stations  F-28l 
and  F-282. 

Comparison  of  data  from  three  depths  at  a  single  station  to 
establish  that  the  indicated  stresses  are  derived  from  a  common  source 
is  not  subject  to  these  limitations  because  ground  ranges  to  all  gages 
at  any  station  are  nearly  identical,  differing  appreciably  only  at 
Station  F-280  for  Shot  9- 

The  general  pattern  of  arrival  times  in  Tables  1.1  and  1.2  is 
roughly  consistent  with  differences  in  ground  range  and  gage  depth. 
However,  arrival  times  are  not  suitably  precise  for  analysis  of  Shot  9 
data  in  terms  of  comparative  propagation  velocities  because  the  obliq¬ 
uity  of  the  array  to  the  air  shock  front  resulted  in  a  very  complex 
system  of  propagation  paths  to  the  gages.  Arrival  times  for  earth 
stresses  from  Shot  10  at  several  stations  show  the  apparent  vertical 
component  of  propagation  velocity  between  the  9-  and  15-Ct  depths  to 
be  about  ’.>00  ft/sec.  This  velocity  component  is  reasonable  Tor  the 
soil  conditions  at  Frenchmen  Flat  and  the  relative  positions  of  che 
gages . 

Rise  times  of  the  earth  stresses  from  'dhots  y  and  10  increase 
directly  with  depth  and  inversely  with  incident  overpressure.  Kxeop- 
bions  to  these  trends  are  usually  in  the  direction  of  greater  rise 
times  and  may  be  a  consequence  of  errors  in  record  interpretation. 

Rise  times  for  the  1-fb  gages  are  short  and  are  probably  characteristic 
of  the  incident  air  shock,  having  been  affected  to  a  negligible  degree 
by  traversing  i  ft  of  soil.  The  increase  in  rise  time  with  depth  is 
consistent  with  transmission  of  stress  through  a  medium  in  which  plastic 
dofnrriiitlnn  occurs.  The  increased  rise  time  as  a  function  of  magnitude 
of  incident  air  overpressure  is  consistent  with,  although  somewhat 
greater  than,  the  rise  time  of  the  air  shock. 

Poativs-phase  duration  in  the  air  shock  wave  for  I)  hot  9  is 
greater  than  in  the  measured  earth  stress.  The re  is  also  a  general 
decrease  of  positive  phase  duration  with  depth.  However,  these  data 
have  the  same  lack  of  cor.gruity  as  do  arrival  times.  Positive-phase 
durations  for  Idiot  10  follow  a  clearer  pattern  of  increase  with  ground 
range  and  decrease  with  depth  except  for  anomalies  in  data  from  the  [>- 
Lt  gages  at  stations  F-28O,  F-216,  and  F-283,  and  in  the  air  overpres¬ 
sure  data  from  Station  F-2ol.  Durations  of  the  precursor  shock  wave 
from  Shot  10  follow  a  pattern  roughly  similar  to  that  for  positive- 
phase  duration  out  to  a  ground  range  of  about  1U0O  ft;  beyond,  it  be¬ 
comes  difficult  to  distinguish  the  precursor  from  the  main  shock  wave. 

Peak  stress-ground  range  curves  for  the  19-ft  gages  at  all  sta¬ 
tions  (Figs.  1.1  and  1.4)  pnrallel  the  trend  established  by  the  data 
from  the  incident  air  shock  wave.  Data  from  the  1-  arid  9-ft  gages  were 
too  erratic  to  define  clearly  any  similar  trends.  Observation  of  earth 
stresses  greater  than  the  incident  air  overpressures  has  not  been 


SECRET  -  RESTRICTED  DATA 


18 

SECRET  -  RESTRICTED  DATA 


*Indetermiaate 


uncommon  in  measurements  of  vertical  earth  stress  at  shallow  depths. 
However,  the  erratic  data  from  the  shallower  gages,  which  contrast 
strongly  with  the  coherent  results  from  those  at  15  ft,  render  the 
high  measured,  stresses  less  dependable  than  the  data  from  the  deep 
gages.  If  the  indicated  peak  streases  for  the  1-ft  gages  at  Stations 
F -216  and  F-282  are  assumed  to  be  low  by  10-20  psi  for  both  Shots  9 
and  10,  a  smooth  curve  through  the  remaining  data  conforms  roughly  to 
that  for  the  air  pressure  data,  as  illustrated  by  the  light  dashed 
curves  of  Figs.  1*3  and  1.4.  A  similar  curve  has  been  drawn  for  the 
5-ft  data  from  Shot  10  by  assuming  that  the  indicated  peak  stresses 
for  that  depth  at  Stations  F-282  and  F-200  are  too  great.  Such  ex¬ 
cessively  high  stresses  may  be  the  result  of  overcompaction  or  exces¬ 
sive  rigidity  of  the  backfill.  Conversely  the  low  stresses  indicated 
by  the  1-ft  data  at  these  stations  and  by  the  5-ft  gage  at  Station 
F-2l6  may  have  been  caused  by  high  compressibility  of  the  Boil  or 
shrinkage  of  the  backfill.  The  minor  scatter  in  the  data  from  15-ft 
gages  suggests  either  that  the  compacted  backfill  was  more  uniformly 
placed  over  them  or  that  its  relatively  small  volume  compared  with 
that  of  the  undisturbed  earth  above  the  gages  obscured  any  irregulari¬ 
ties  introduced  by  the  backfill. 

Data  on  arrival  times  and  peak  stresses  suggest  that  the  air 
shack  Incident  at  the  station  is  the  source  of  the  observed  earth 
stresses.  A  stronger  indication  that  air  overpressure  is  the  major 
source  of  the  earth  stress  is  given  by  superposition,  with  coincident 
arrival  times,  of  the  air  overpressure  curve  from  Station  F-28l  and 
the  earth  stress  curves  for  the  1-  and  5-ft  gages  at  Station  F-282  for 
Shot  10  (Fig.  1.5).  The  pattern  registered  by  the  earth  stress  gages 
is  very  similar  to  that  registered  by  the  air  pressure  gage  except  for 
discrepancies  in  magnitude.  Stress-time  curves  of  Figs.  A.l  through 
A. 5  also  reveal  similarity  between  curves  for  the  gages  of  each  group, 
and  characteristic  features  of  the  earth  stress  curves  Eire  also  found 
in  air  overpressure  curves  from  BtationB  too  remote  for  direct  compari¬ 
son  . 

The  1-ft  gages  at  Stations  F-282  and  F-200  indicated  lower  peak 
stresses  than  the  5-ft  gages.  This  condition  is  a  reversal  of  the 
trend  at  the  other  three  stations  and  is  contradicted  by  data  from  the 
15-ft  gages  at  all  five  stations.  Rise  times,  arrival  times,  and  posi¬ 
tive-phase  durations  support  the  data  as  tabulated .  Consequently  this 
reversal  in  stress  magnitudes  must  be  considered  real  and  probably  the 
result  of  discrepancies  in  gage  and  backfill  placement. 

No  essential  difference  in  performance  between  older  find  newer 
gages  is  evident  from  the  stress-time  curves  (Figs.  A. 2,  A. 4,  and  A. 5)* 
The  5-ft  gage  (old)  did  register  considerably  higher  stresses  than  the 
1-ft  gage  (new)  at  Station  F-282,  but  a  similar  situation  was  observed 
at  Station  F-200  where  both  gage3  were  new  ones  (Figs.  1.3  and  1.4). 
Data  from  the  15-ft  gage  (old)  at  Station  F-283  deviated  less  from  the 
curves  of  peak  stress  vs  ground  range  them  did  data  for  new  gages  at 
two  other  stations. 

There  seems  to  be  no  clear-cut  general  evidence  of  ground-trans¬ 
mitted  stress  from  remote  ground  incidence  of  the  air  shock.  It  is 
possible  that  the  minor  signal  at  about  1.5  sec  in  the  record  from  the 
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Fig.  1.5  Comparison  of  Air  Overpressure  mid  Verticnl  Earth  Stress  - 
Shot  10 


l'p-ft  gage  at  Station  F-280  for  Shot  9  (Fig.  A.l)  may  represent  re¬ 
fracted.  or  reflected  earth  stress  originating  from  incidence  of  the 
air  shock  in  the  vicinity  of  ground  zero.  This  signal  is  repeated  at 
the  shallower  gages  with  a  finite  time  lag  and  decreased  amplitude, 
suggesting  transmission  upward.  A  similar  phenomenon  appears  in  the 
records  from  Station  F-283  for  Shot  9  (Fig.  A.2),  but  no  correspond¬ 
ingly  distinct  evidence  appears  at  other  stations  for  this  shot.  Data 
from  Shot  10  give  no  evidence  of  similar  ground-transmitted  signals, 
perhaps  because  they  are  obscured  bv  perturbations  introduced  by  the 
precursor  shock.  The  secondary  positive  signals  which  appear  about 
Go  sec  after  main  shock  arrivals  from  Shot  10  at  Station  F-280  (Fig. 

A. 3)  and  at  progressively  greater  intervals  after  main  shock  arrival 
at  more  remote  stations  are  the  result  of  a  secondary  shock  in  the  air 
overpressure.  The  increasing  time  interval  between  the  main  and  sec¬ 
ondary  events  at  increasing  ground  ranges  is  not  consistent  with  arriv¬ 
al  of  n  secondary  event  through  the  earth  from  the  vicinity  of  ground 
zero. 

The  effect  of  gage  depth  on  observed  stresses  has  been  investi¬ 
gated  by  plotting  stress  magnitude  against  depth.  Data  for  this  study 
were  read  at  points  judged  to  represent  common  events  on  the  stress¬ 
time  curves  for  a  single  station.  Several  comparable  sets  of  data, 
including  peak  stresses,  were  derived  from  each  group  of  gages.  Such 
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events  did  not  occur  simultaneously  at  all  depths  hut  lagged  at  the 
deeper  gages  by  intervals  consistent  with  the  corresponding  lags  in 
arrival  time.  Stress  varied  as  an  inverse  function  of  depth,  but  the 
variation  was  not  linear.  Tests  were  made  to  determine  whether  attenu¬ 
ation  could  be  better  expressed  by  an  equation  of  the  form 


P  =  T1  exp  -  p  (d  -  dx)  (1.1) 

or 

P  =  Px  (d  /  d^-®,  (1.2) 


where  P  and  are  stresses  which  exist  at  depths  d  and  d]_  and  £  and  a 
are  constants . 

Semilogarithmic  plots  of  all  sets  of  comparable  data  fox*  each 
station  were  made,  and  a  linear  relation  for  each  set  was  sought.  Data 
for  all  three  depths  fell  very  nearly  on  a  straight  line  in  some  in¬ 
stances,  and  in  general  they  could  be  roughly  approximated  by  a  straight 
line.  Peak  stress  data  usually  approximated  linearity  about  as  well  as 
did  data  representing  other  events  in  the  stress-time  curves.  There¬ 
fore  they  were  used  as  a  basis  for  final  analysis  of  data  from  each 
station.  However,  the  results  may  be  considered  to  represent  the  stress- 
depth  relation  of  the  whole  stress  transient  rather  than  of  peak  stresses 
only. 

Semilog  and  log-log  plots  of  peak  stress  versus  depth  for  Shot  9 
are  presented  in  Figs.  1.6  and  1.7-  Corresponding  plots  for  Shot  10 
arc  presented  as  Figs.  1.0  and  1.9*  The  coherence  and  negligible  scat¬ 
ter  of  data  from  the  15-ft  gage3  for  both  shots  Indicated  that  they 
were  the  most  trustworthy  and  should  be  given  maximum  weight.  Con¬ 
versely  data  from  some  of  the  5-ft  gages,  particularly  at  Station  F-282 
and  F-200,  were  ignored  because  of  their  anomalous  nature.  Data  from 
Shot  9  for  the  5-ft  gages  at  Station  F-2l6  and  from  Shot  10  for  the  5~ 
ft  gages  at  Station  F-280  were  similarly  ignored,  the  former  being  ex¬ 
cessively  low,  the  latter  anomalously  high.  Data  from  the  1-ft  gages 
were  given  full  weight  although  this  procedure  was  sometimes  arbitrary 
and  resulted  in  inconsistently  flat  slopes.  The  scarcity  of  data  and 
the  large  proportion  of  Intermediate  data  points  that  were  ignored  be¬ 
cause  of  anomalies  left  little  choice  in  the  lines  drawn  to  represent 
the  stress-depth  relation  since  only  two  points  remained. 

.Ml  data  from  several  stations  fell  on,  or  nearly  on,  a  straight 
line  on  one  or  the  other  of  the  graphs .  Data  from  Station  F-203  fit  a 
straight  line  on  the  log-log  plots  for  both  Shots  9  and  10.  Data  from 
Station  F-280  for  Shot  9  and  from  Station  F-216  for  Shot  10  fell  on  a 
straight  line  on  the  semilog  plot.  Several  other  sets  of  data  are 
reasonably  well  represented  by  straight  lines.  Analytic  plots  of  data 
from  Station  F-202  for  Shot  9  substantiate  evidence  from  Fig.  1.2  that 
they  were  too  irregular  to  yield  a  logical  relation  between  stress  and 
depth.  Data  from  this  station  for  Shot  10  are  more  reasonable,  but 
again  they  are  the  most  irregular  of  the  five  sets  of  information  avail¬ 
able  . 


22 

SECRET  -  RESTRICTED  DATA 


Cape  Depth  (ft) 


Vertical  Earth  Stress 
Attenuation  -  Shot  9 


,8  Vertical  Earth  Stress 
Attenuation  -  Shot  10 
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TABLE  1.3  -  Attenuation  Factors 


Station 

_ fKf  t  ) _ 

1 

Shot  9 

Shot  10 

Shot  10 

Shot  9 

Shot  10 

number 

F-280 

19 

140 

0.058 

0.39 

0.30 

F-216 

17 

72 

0 .092 

0.0h6 

0.26 

0.l4 

F-262 

16 

33 

0.016 

O.O63 

0.09 

0.30 

F-283 

19 

0.077 

0.071 

0.40 

O.36 

F-200 

_ 

12 

8 

0.098 

O.083 

0.49 

0.45 

Mean 

O.O76 

o.o64 

0-39 

0.35 

Attenuation  factors  for  Shot  9  at  Station  F-282  and  for  Shot  10  at  Sta¬ 
tion  F-2l6  were  unusually  low  and  have  been  neglected  in  computing  the 
mean  values  in  Table  1-3*  There  is  a  definite  suggestion  of  increasing 
values  of  p  and  a  with  decreasing  peak  incident  overpressure  (increas¬ 
ing  ground  range)  for  Shot  10;  no  trend  of  this  sort  is  evident  for 
Shot  y.  This  apparent  inconsistency  is,  however,  a  consequence  of  the 
relative  spread  of  peak  overpressure  represented*  The  peak  overpres¬ 
sures  incident  at  the  last  two  stations  for  Shot  10  are  commensurate 
with  the  full  range  of  overpressures  for  Shot  9>  and  the  attenuation 
factors  for  Shot  10  at  these  two  stations  are  roughly  comparable  with 
those  for  Shot  9  at  all  stations .  Scatter  of  the  basic  data  and  the 
somewhut  arbritary  procedures  adopted  in  their  analysis  preclude  use  of 
the  resulting  factors  for  more  than  rough  estimates  based  on  average 
values  or  indication  of  trends.  The  analysis  shows  that  attenuation  of 
vertical  earth  stress  as  a  function  of  depth  may  be  expressed  equally 
well  by  the  equation 

P  =  P]_  exp  [-  0.07  (d  -  dx)  ]  , 

or  by  the  equation 

P  =  PL  (d  /  d^"0*37 

where  depths  are  measured  in  feet.  The  factor  p  is  probably  accurate 
within  +  20  per  cent  and  the  exponent  a  within  +  25  per  cent. 

It  is  comforting  to  note  the  close  agreement  between  the  vuiues 
of  p  and  a  as  calculated  from  data  fur  Shots  9  and  10.  It  is  somewhat 
disturbing,  however,  that  the  fit  of  the  data  is  no  better  for  one 
form  of  analysis  than  for  the  other.  Results  of  the  two  forms  of  anal¬ 
ysis  as  shown  in  Figs.  1.6  through  1.9  certainly  do  not  provide  a  basis 
for  regarding  one  type  of  attenuation  more  favorably  than  the  other. 
There  are  intuitive  grounds  for  preferring  the  logarithmic  type  of 
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attenuation,  but  In  the  absence  of  stronger  experimental  evidence  for 
differentiation  between  the  logarithmic  and  inverse  power  types  of  at¬ 
tenuation,  intuition  hardly  seems  sufficient  grounds  for  tipping  the 
balance. 

Integration  of  the  stress -time  curves  gives  curves  of  impulse  in¬ 
tensity  or  impulse  per  square  inch  as  a  function  of  time  for  each  gage 
at  the  five  stations.  Families  of  curves  for  each  station  are  pre¬ 
sented  in  Appendix  A  as  Figs.  A. 6  through  A. 9*  Integration  was  per¬ 
formed  by  summation  of  trapezoidal  areas  derived  from  original  record 
data  and  entered  on  punched  cards .  IBM  machines  were  used  for  the 
mathematical  operations .  No  correction  was  made  for  zero  drift  or 
residual  stresses  which  were  evident  in  some  of  the  recorded  data. 
Errors  from  these  sources  do  not  become  excessive  during  the  first 
portion  of  tiie  integrated  curve  and  probably  do  not  affect  the  maxima 
by  more  than  10  per  cent.  Beyond  the  maxima  the  integrated  curves  are 
strongly  affected  by  the  cumulative  effect  of  residual  stresses,  but 
these  latter  portions  of  the  curves  are  of  little  interest. 

The  curves  of  impulse  intensity  versus  time  and  the  peak  impulse 
intensities  listed  in  Tables  1.1  and  1.2  are  presented  for  information 
only.  No  special  significance  is  attached  to  them.  Variation  of  im¬ 
pulse  with  depth  is  analogous  to  that  of  earth  stress  but  more  pro¬ 
nounced  because  both  stress  magnitude  and  positive-phase  duration  de¬ 
crease  with  depth.  Curves  for  Shot  10  show  the  effect  of  the  precursor 
shock  on  impulse  intensity.  The  precursor  makes  an  appreciable  con¬ 
tribution  because  of  its  duration  even  though  it  involves  relatively 
low  stresses.  The  rise  in  the  precursor  portion  of  the  curve  is  not 
steep,  but  its  contribution  to  the  maximum  impulse  is  equivalent  to  a 
high  peak  stress  of  a  nonprecursor  shock  wave  since  such  a  peak  would 
have  a  short  duration  compared  with  that  of  a  well  developed  precursor 
shock. 


i.6  CONCLUSIONS  AND  RECOMMENDATIONS 


1.  Observed  earth  stresses  were  induced  by  incidence  of  air 
shock  at  the  ground  surface  in  the  vicinity  of  the  gages . 

2.  The  vertical  component  of  earth  stress  so  induced  decreases 
with  depth.  The  attenuation  may  be  expressed  equally  well 
by  the  equations 

P  =  ?1  exp  [-  0.07  (d  -  dx)] 

or 

p  =  p1  (d  /  d1)'0,37  , 

where  P  and  P^  are  the  peak  stresses  in  psi  at  depths  d  and 
d]_  in  ft.  The  constants  in  both  equations  are  probably  good 
within  +  25  per  cent. 

3.  Stress  magnitudes  of  the  entire  stress-time  sequence  are 
attenuated  in  essentially  the  same  manner  as  peak  stresses. 
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4.  Data  from  Shot  10  suggest  that  rate  of  attenuation  increases 
with  decreasing  incident  overpressure.  This  increase  appeal's 
to  be  roughly  logarithmic,  hut  the  data  are  not'  adequate  to 

ccnd.usi.vc  • 

5.  Rise  time  of  earth  stress  increases  and  positive -phase  dura¬ 
tion  decreases  with  depth.  Arrival  times  lag  with  increasing 
depth . 

6.  Peak  impulse  intensities  vary  with  depth  analogously  to,  but 
more  rapidly  than,  the  stresses. 

Unsolved  instrumentation  problems  and  incomplete  knowledge  of  the 
phenomena  of  transmission  of  air-shock-induced  earth  stress  cause  the 
results  of  these  measurements  to  be  subject  to  controversy.  The  situa¬ 
tion  will  not  be  clarified  by  more  tests  of  the  same  sort  because  these 
uncertainties  are  fundamentally  the  result  of  instrumentation  problems. 
The  following  recommendations  are  presented  as  a  basis  for  removing 
these  uncertainties : 

1.  Effects  on  measured  stresses  of  instrument  characteristics 
and  of  materials  and  methods  of  gage  placement  should  he 
studied.  Such  studies  should  he  made  initially  in  soils 
typical  of  test  sites  but  ultimately  should  include  soils 
typical  of  target  areas.  The  goal  of  these  tests  would  be 
firm  knowledge  of  the  reliability  of  indicated  stresses  in 
terms  of  stresses  in  the  undisturbed  medium.  Static  loads 
and  dynamic  loads  induced  by  high-explosives  charges  of 
moderate  size  would  be  used  for  these  tests.  Cognizance 
would  have  to  be  taken  of  possible  adverse  effects  of  the 
short  time  factor  inherent  in  small  charges  since  the  goal 
of  the  instrumentation  is  use  under  conditions  produced  by 
very  large  charges. 

2.  Full-scale  free-field  measurements  similar  to  those  of  Pro¬ 
ject  1.4a  should  he  performed  at  Nevada  Proving  Ground  after 
instrumentation  problems  have  reached  a  reasonable  solution. 
The  full-scale  experiment  should  include  earth  stress  and 
earth  strain  measurements  at  several  depths,  including  the 
maximum  feasible  for  gage  placement.  Incident  overpressure 
levels  should  range  from  300  to  30  psi.  Huch  tests  should 
include  means  for  observing  the  effects  of  degree  of  com¬ 
paction  of  the  soil  on  possible  attenuation  of  stresses. 

The  ultimate  purpose  of  free-field  measurements  of  earth  stress 
attenuation  is  corrslation  between  free-field  phenomena  and  structural 
response  or  damage  to  underground  facilities.  Such  tests  must  conse¬ 
quently  be  closely  related  to  experiments  of  the  type  performed  by 
Project  3.8  of  UPSHOT -KNOTHOLE,  and  the  recommendations  contained  in 
the  final  report  on  Project  3-8  tests  of  buried  structures  arc  con¬ 
sidered  desirable. 

Theoretical  studies  aimed  specifically  at  better  understanding  of 
the  transmission  and  dissipation  of  steep-fronted  transient  stresses  in 
soils  would  aid  materially  in  interpreting  free-field .data  from  stress 
measurements . 
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PART  IT 


DERIVATION  OF  EARTH  STRESS-STRAIN  SYSTEMS  AND  GROUND-MOTION  STUDIES 


2.1  THEORETICAL  BACKGROUND 


Damage  to  an  underground  structure  as  a  result  ol  an  explosion 
will  be  directly  related  to  the  transient  force  field  or  one  of  its  de¬ 
finitive  parameters ,  and  will  result  from  the  response  of  the  structure 
or  its  elements  to  one  or  more  of  these  phenomena.  Consequently  the 
problem  of  predicting  damage  to  underground  structures  is  more  complex 
than  simply  defining  the  force  field  produced  in  the  earth  by  an  ex¬ 
plosion.  In  fact,  it  is  resolved  Into  problems  of  (.1)  defining  the 
force  field  by  free-field  parameters,  (2)  interaction  of  the  structure 
with  the  force  field,  and  (3)  determining  consequent  loading  and  struc¬ 
tural  response.  Structural  tests  of  the  type  undertaken  by  Project  3-^ 
of  UPSHOT- KNOTHOLE^/  will  give  answers  to  the  third  problem  and  may 
help  clarify  the  second.  Specific  cases  involving  a  particular  struc¬ 
ture  at  a  known  distance  from  an  explosion  of  known  magnitude  can  bo 
treated  by  direct  measurement  of  structural  response.  However,  this 
approach  must  necessarily  be  of  limited  scope,  and  the  general  approach 
would  be  far  more  useful  if  free-field  phenomena  were  understood  and 
could  be  correlated  with  structural  response  and  damage. 

Definition  of  the  transient  force  field  in  the  absence  of  a  struc¬ 
ture  requires  knowledge  of  the  stress  and  strain  tensors  and  the  accel¬ 
eration,  velocity,  or  displacement  of  soil  particles  at  various  distances 
from  the  explosion  and  at  various  instants  of  time  following  detonation. 
Knowledge  of  all  three  parameters  --  stress,  strain,  and  motion  --  is 
required  because  structures  of  different  types  will  be  affected  more 
strongly  by  one  or  the  other.  Initially  the  feasibility  of  measuring 
these  parameters  in  the  free  field  must  be  established.  Ultimately  such 
measurements  as  are  feasible  must  be  made  in  the  free  earth  and  in  the 
earth  adjacent  to  structures  and  the  results  correlated  with  measured 
loading  and  response  of  structural  elements. 

2.1.1  Earth  Stress  Systems 

Derivation  of  principal  stresses  and  strains  and  of  stress  and 
strain  tensors  is  treated  in  considerable  detail  in  textbooks  on  theory 
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of  elasticity  and  stress  analysis.®/  These  sources  are  usually  con¬ 
cerned.  with  static  conditions ,  uus  tne  meunocis  ox  analysis  can  "be  ap¬ 
plied  to  instantaneous  values  of  transient  phenomena. 

The  stress  tensor  at  a  point  in  a  medium  is  defined  by  the  mag¬ 
nitude  and  directions  of  the  normal  and  shear  stresses*  on  any  group 
of  three  orthogonal  plenes  at  the  point.  This  definition  requires 
measurement  of  at  least  six  parameters,  three  normal  stresses  and  three 
shear  stresses.  An  equally  valid  and  often  more  useful  definition  of 
the  stress  tensor  is  in  terns  of  magnitudes  and  directions  of  the  prin¬ 
cipal  stresses;  these  are  the  normal  stresses  on  the  three  principal 
planes,  which  are  mutually  orthogonal  and  so  oriented  that  all  shear 
stresses  on  them  vanish.  The  normal  and  shear  s  tresses  on  any  plane 
through  the  common  point  of  the  principal  planes  can  be  derived  from 
the  principal  stresses. 

Techniques  and  procedures  for  measurements  from  which  the  stress 
tensor -may  be  derived  have  been  developed  by  the  Corps  of  Engineers,  US 
Army, 2/  in  the  course  of  studies  of  stress  distribution  in  earth-fill 
structures  and  externally  loaded  runway  subgrades .  In  these  studies 
loading  was  static,  and  either  plane  or  radially  symmetric  stress  could 
be  assumed. 


Measuring  transient  earth  stresses  introduces  several  problems 
not  encountered  in  static  measurements.  Time  dependence  of  loading  and 
its  consequences,  such  as  altered  plastic  behavior  of  the  soil  and  in¬ 
ertial  effects,  are  the  most  evident  differences.  Instrumentation  prob¬ 
lems  differ  appreciably,  although  the  relative  difficulty  of  static  and 
dynamic  measurements  may  be  debatable.  Instruments  suitable  for  meas¬ 
uring  shear  stresses  in  soils  under  static  load  do  not  appear  to  be 
practical  for  use  under  transient  loads,  and  shear  stresses  must  conse¬ 
quently  be  derived  from  measured  normal  stresses.  Assumption  of  either 
a  plane  or  radially  symmetrical  stress  system  is  less  likely  to  be 
permissible  under  dynamic  them  under  static  loading.  Consequently  meas¬ 
urements  under  dynamic  conditions  should  be  designed  to  furnish  suffi¬ 
cient  component  data  to  permit  analysis  as  a  three-dimensional  stress 
system. 


Stress  tensors  derived  from  earth  stress  measurements  during 
dynamic  loading  are  customarily  for  dynamic  stresses  only  and  do  not 
define  the  total  stress  condition  because  the  tensor  for  static  earth 
stress  is  ignored  in  both  observation  and  analysis.  This  situation 
results  principally  from  differences  in  instrumentation  for  the  two 
classes  of  measurements.  The  static  stress  tensor  may  be  determined 
from  independent  static  measurements  either  before  or  after  dynamic 
loading.  The  complete  stress  condition  is  then  defined  by  the  sum  of 
the  static  and  dynamic  tensors. 

Numerous  measurements  of  earth  pressure  resulting  from  under¬ 
ground  explosions  have  been  made  by  means  of  oil-bag  pressure  gages. 12/ 
The  value  of  these  measurements  is  uncertain  because  of  their  doubtful 


*  Normal  stress  on  a  plane  is  the  force  intensity  or  force  per  unit 
area  acting  perpendicular  to  the  plane.  The  corresponding  shear 
stress  is  the  force  intensity  acting  parallel  to  the  plane. 
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relation  to  stresses  and  their  critical  dependence  on  such  factors  as 
tne  column  of  fluid  in  the  gage  hole.  They  proved  useful,  however,  in 
providing  a  basis  for  correlation  between  phenomena  produced  by  explo¬ 
sions  of  various  magnitudes  in  several  types  of  soil.i= 

The  first  tests  of  instrumentation  to  measure  transient  earth 
stress  were  made  in  conjunction  with  the  underground  high-explosive 
tests  of  Project  1  (9)  of  JANGLE. iS/  At  the  invitation  of  Stanford 
Research  Institute,  Sandia  Corporation  installed  instrumentation  to 
measure  four  components  of  earth  stress  on  Shots  HE-1  and  HE -3 ,i3/ 
Carlson-Wiancko  earth  stress  gages  were  placed  5  ft  deep,  oriented  to 
respond  to  normal  stresses  on  the  vertical,  radial,  and  two  45°  vertical 
radial  planes*  at  a  specific  ground  range.  The  preliminary  nature  of 
these  tests  and  the  limited  time  and  number  of  recording  channels  avail¬ 
able  prompted  assumption  of  a  radially  symmetrical  stress  system.  Re¬ 
sults  indicated  that  measuring  techniques  for  static  earth  stress 
components  could  be  adapted  to  dynamic  stress  measurements.  Certain 
alterations  of  instrument  deslgn2/  were  found  necessary,  however,  to 
eliminate  undesirable  mechanical  resonance  (ringing)  of  the  gages. 

More  recently  Stanford  Research  Institute  has  measured  earth 
stress  components  as  a  part  of  Project  Mole. 12/  Unidirectional  stress 
gages  were  oriented  within  the  earth  to  respond  to  normal  stresses  on 
the  vertical,  radial,  and  tangential  planes  as  well  as  on  several  of 
the  45°  bisecting  planes  of  the  orthogonal  group.  Check  tests  included 
duplication  of  gages  at  several  orientations  as  well  as  variation  of 
depth  and  placement  conditions.  Stress  component  data  were  insuffi¬ 
cient  for  complete  three-dimensional  analysis,  but  considerable  infor¬ 
mation  pertinent  to  transient  earth  stress  problems  was  obtained. 

2.1.2  Earth  Strain  Systems 

Knowledge  of  the  strain  tensor  at  a  point  in  the  earth  is  im¬ 
portant  to  understanding  underground  damage  phenomena.  Correlation 
of  stresses  and  strains  as  a  function  of  time  gives  needed  information 
concerning  the  elastoplastic  behavior  of  soil  under  transient  loading. 
Moreover,  certain  categories  of  structures  are  more  sensitive  to  strain 
phenomena  in  the  earth  than  to  stress  loading. 

Static  earth  strain  components  have  been  measured,  as  a  part  of 
stress -distribution  studies  of  loaded  runway  subgrades. 2/  Strain  gages 
developed  at  the  Ohio  River  Division  Laboratory  (ORDL)  of  the  Corps  of 
Engineers,  US  Army,  have  proved  satisfactory  for  measurement  of  static 
earth  strains.  However,  results  obtained  by  Stanford  Research  Insti¬ 
tute  with  ORDL  strain  gages  subjected  to  transient  loading  during  Phase 
I  of  Project  Mole  were  not  wholly  satisfactory  although  more  recent  in¬ 
formation  has  tempered  initial  adverse  evaluation. iQ' 


*  The  direction  of  a  plane  is  defined  as  the  direction  of  the  normal 
to  it  throughout  this  report,  and  consequently  corresponds  to  the 
direction  of  the  normal  stress  which  acts  on  it. 
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Sandia  Laboratory  developed  an  earth  strain  gage,  similar  in 
principle  to  the  OKDL  gage,  for  Project  I9.U1  of  TUMBLER-SNAPPER. 3/ 

It  was  designed  to  have  greater  sensitivity,  smaller  maximum  range, 
and  density  more  closely  approximating  that  of  the  soil)  characteris¬ 
tics  considered  essential  for  proper  response  to  transient  earth  strains. 
Tests  of  these  Sandia  gages  during  TUMBLER -SNAPPER  showed  them  to  be 
suitable  in  principle  but  indicated  that  transducers  of  greater  elec¬ 
trical  sturdiness  were  needed  to  survive  the  zero-time  electrical  tran¬ 
sients  associated  with  nuclear  detonations. 

2.1.3  Earth  Accelerations 


Measurement  of  one  or  more  mutually  orthogonal  components  of 
earth  acceleration  furnishes  information  from  which  motion  of  the  earth 
and  of  the  underground  portions  of  some  structures  may  be  derived  for 
correlation  with  damage.  Integration  of  the  acceleration-time  data 
provides  velocity-time  data  which  may,  in  turn,  be  integrated  to  give 
displacement-time  infonnation.  Accelerometers  axe  the  most  suitable 
instruments  for  measurement  of  earth  motion  produced  by  nuclear  explo¬ 
sions.  Magnitudes  and  durations  of  velocities  and  displacements  are 
so  great  that  development  of  velocity  or  displacement  gages  has  not 
been  practicable. 

Earth  acceleration  components  resulting  from  air  and  under¬ 
ground  nuclear  and  high-explosives  bursts  have  been  measured  by  numer¬ 
ous  agencies. ili/  Instrumentation  techniques  are  well  established,  and 
data-handling  procedures  are  comparatively  simple.  Patterns  of  recorded 
accelerations  have  not  been  generally  simple  or  consistent.  In  fact, 
about  the  only  consistently  dependable  feature  of  these  measurements 
has  been  the  strong  initial  downward  acceleration  that  characterizes 
arrival  of  the  air-shock-induced  Bignal  at  the  vertical-component  ac¬ 
celerometer.  Reasonable  displacement  information,  derived  from  double 
integration  of  the  acceleration-time  records,  has  been  reported  in  some 
instances.  Displacement  data  so  derived  are,  of  course,  subject  to  the 
errors  inherent  to  integration  and  correction  for  variation  of  record 
zero,  but  data  which  permit  such  manipulation  without  obviously  gross 
errors  produce  reasonable  information. 

2.2  PLAN  OF  THE  EXPERIMENT 


Loading  of  the  earth  by  air  shock  from  nuclear  air  bursts  provides 
excellent  conditions  for  testing  and  verifying  both  instrumentation  and 
analytical  procedures  for  determining  earth  stress  tensors.  Such  .load¬ 
ing  induces  a  stress  tensor  characterized  by  a  major  principal  stress 
directed  very  nearly  vertically  and  considerably  greater  than  the  other 
two  principal  stresses.  Magnitudes  and  durations  of  load-induced  phe¬ 
nomena  would  be  sufficient  to  minimize  the  effect  of  the  finite  size 
of  an  instrument  array  on  the  assumption  that  all  measured  stresses 
represent  data  at  a  single  point  within  the  earth.  Neither  of  these 
conditions  was  realized  in  tests  prior  to  UPSHOT-KNOTHOLE. 

Project  1.4b  was  conceived,  therefore,  as  a  test  of  instrumen¬ 
tation  and  analysis  for  derivation  of  induced  earth  stress  tensors. 
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mination  of  the  tensor  and  were  to  be  supplemented  by  teste  of  improved 
Sandia  earth  strain  gages  and  routine  measurements  of  acceleration  com¬ 
ponents.  Measurements  were  planned  for  three  shots  of  UPSHOT-KNOTHOLE, 
and  all  instrumentation  was  to  be  placed  at  a  depth  of  5  ft  except  that 
for  ground -baffle  measurements  of  air  overpressure  at  the  surface. 

Instrumentation  for  Shot  1  provided  for  observation  at  the  an¬ 
ticipated  20-psi  air  overpressure  ground  range  of  six  stress  components, 
the  corresponding  six  strain  components,  three  acceleration  components, 
and  incident  air  overpressures  (Fig.  2.1).  Duplicate  instrumentation 
was  included  to  provide  a  reliability  chock  on  the  measurements  of  two 
stress  and  one  strain  component.  Tests  planned  for  Shots  9  and  10  at 
Frenciiman  Flat  were  more  comprehensive  than  for  Shot  1  and  included 
instrumentation  for  measuring  nine  stress,  six  strain,  and  three  ac¬ 
celeration  components  in  the  earth  as  well  as  incident  air  overpressure 
(Fig.  2.2).  Multiple  instrumentation  was  planned  for  two  stress  and 
one  strain  component  for  reliability  checks.  Identity  of  nominal  ground 
zeros  for  Shots  9  and  10  permitted  use  of  a  single  instrument  array  for 
both  tests,  but  differences  in  burst  height  and  yield  required  the  use 
of  set  ranges  fur  anticipated  pc  ale  incident  overpressure  of  20  psi  l'or 
Shot  9  and  30  psi  for  Shot  10.  Limited  recording  facilities  necessi¬ 
tated  measurement  of  only  nine  stress  components,  including  one  dupli¬ 
cation,  for  Shot  9*  However,  all  instrumentation  for  this  project  was 
operated  during  Shot  10. 


2 . 3  INSTRUMENTATION 

Instrumentation  used  for  Project  1.4b  measurements  was  generally 
similar  to  that  used  during  previous  tests  at  Nevada  Proving  Ground. 

The  Carlson-Wiancko  earth  stress  gages  and  Sandia  earth  strain  gages 
were  modified  as  described  below.  Wiancko  accelerometers  and  Wianeko 
air  pressure  gages  were  of  standard  design.  Each  end  instrument  was 
connected  through  a  buried  4 -conductor  shielded  .cable  to  an  Qualifier 
and  recorder  in  a  remote  subterranean  shelter .9/  A1'  ;nd  instruments 

were  calibrated  on  assigned  recording  channels,  and  recorder  set  ranges 
and  calibration  steps  were  adjusted  prior  to  placement  of  the  earth 
gages  in  the  ground.  Circuits  of  the  gages  at  Station  F-201  were  checked 
between  Shots  9  and  10,  but  recalibration  for  Shot  10  was  not  feasible 
except  for  the  ground-baffle  gages. 


2.3.1  End  Instruments 

Modification  of  the  Carlson-Wiancko  earth  stress  gage^/  was  a 
purely  geometric  one  to  improve  its  response  aspect  for  earth  stress. 
Prior  to  modification  it  had  been  proven  an  electrically  and  mechani¬ 
cally  reliable  instrument  for  dynamic  measurements.  However,  the  ap¬ 
parent  compressibility  of  the  unmodified  gage  was  appreciably  lower 
than  that  of  its  load-responsive  flat  cylindrical  portion  because  of 
the  transducer  housing  (Fig.  2.3  left).  Studies  of  the  effect  of  gage 
compressibility  on  the  distribution  of  static  stress  in  soil  have  shown 
that  if  this  compressibility  is  low  compared  with  that  of  the  soil, 
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Fig.  2.2  Instrumentation,  Frenchman  Flat 
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Pig.  2.3  Carlson-Wiancko  Earth  Stress  Gages;  Unmodified 
(left),  Modified  (right) 


concentration  of  stress  occurs  at  the  gage.  The  relative  significance 
of  the  ratio  of  instrument  and  soil  compressibilities  for  transient 
dynamic  versus  static  soil  loads  is  not  fully  understood.  It  is  prob¬ 
able  that  higher  stress-strain  ratios  of  the  soil  under  dynamic  loads 
brings  the  ratio  of  gage  and  soil  compressibilities  closer  to  the  ideal. 
Nevertheless  it  was  considered  advisable  to  improve  gage  geometry  as 
a  means  of  improving  stress  response.  Reorientation  of  the  transducer 
element  at  the  expense  of  a  slight  loss  in  sensitivity  permitted  re¬ 
duction  of  the  height  of  the  transducer  housing  from  4.5  to  0.75  in. 

(Fig.  2.3,  right).  These  gages  were  redesigned  with  the  cooperation 
of  Dr.  Roy  W.  Carlson  and  the  Wiancko  Engineering  Co. 

Earth  strain  gages  for  Project  1.4b  were  similar  to  those  in¬ 
itially  tested  in  the  earth  mound  during  TUMBLER-SNAPPER. 3/  Electrical 
weakness  of  the  prototype  gage  was  overcome  by  use  of  better  transducer 
elements.  Improved  mechanical  characteristics  resulted  from  more  rigid 
end  plates  and  a  longer  sliding  surface  of  reduced  area.  The  gage  length 
was  increased  to  5  in.  to  alleviate  the  crowded  conditions  in  the  trans¬ 
ducer  connector  portion  of  the  prototype  gages.  A  better  density  match 
with  the  adjacent  soil  was  obtained  by  use  of  different  construction 
materials.  The  over-all  density  of  the  improved  gages  was  about  100 
lb/ft3,  approximately  the  density  of  the  compacted  Nevada  Proving  Ground 
soil  in  the  backfill.  . 

Accelerometers— iuii?/  were  mounted  in  a  case  fabricated  from  a 
standard  10-in.  bull  plug.*  The  cover  plate  of  the  case  provided  means 


*  A  bull  plug  is  a  standard  pipe  fitting  for  terminating  a  large-diam¬ 
eter  pipeline. 
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for  mounting  three  accelerometers  with  'response  axes  mutually  orthogonal 
(Pig.  2.4).  The  assembly  was  designed  to  have  a  density  approximately 
equal  to  that  of  the  surrounding  soil.  . 

Duplicate  air  pressure  gagesi^ulS/  were  placed  in  a  concrete 
block  ground  baffle  adjacent  to  the  earth  gages  to  indicate  the  inci¬ 
dent  overpressure  at  the  ground  surface. 

2.3.2  Pl.acement  of  End  Instruments 

Instrumentation  for  Project  1.4b  was  installed  at  Station  S-285 
(Area  T-3,  1450  ft  from  ground  zero)  for  Shot  1  and  at  Station  F-28l 
(Frenchman  Flat,  1200  ft  from  nominal  ground  zero)  for  Shots  9  and  3.0. 
All  instrumentation  except  the  air  overpressure  gages  in  ground  baffles 
was  offset  normally  from  the  blast  line  at  the  specified  stations. 

Earth  stress  and  strain  gages  were  placed  in  trenches  5  ft  deep  having 
a  bottom  width  of  about  2  ft  and  a  width  of  about  7  feet  at  the  surface. 
Stress  gages  were  placed  at  3-ft  intervals,  and  strain  gages  were  ap¬ 
proximately  midway  between.  Orientation  was  controlled  within  +3  de¬ 
grees,  and  gage  elevations  were  held  within  +0.1  ft  of  the  prescribed 
direction  and  depth.  Backfill  materials  and  controls  were  similar  to 
those  described  in  Part  I,  Section  1.3*  Details  of  gage  layout  and 
placement  procedures  are  discussed  in  Appendix  B. 

The  accelerometer  assemblies  were  placed  in  holes  5  ft  deep 
and  12  in.  in  diameter.  Accelerometers  were  oriented  with  respect  to 
nominal  ground  zero  by  means  of  the  sight-tube  jig  described  in  Appen¬ 
dix  B.  The  assembly  was  bonded  to  the  soil,  after  orientation,  by  a 
Calseal  grout  and  the  hole  backfilled  with  tamped  soil. 

2.4  RESULTS 

Raw  data  from  Project  1.4b  comprised  oscillograph  and  magnetic 
tape  x'ecords  of  the  various  gage  output  signals  as  a  function  of  time. 
These  records  were  read  and  translated  into  stresses,  strains,  and  ac¬ 
celerations  by  automatic  computers  and  were  furnished  for  analysis  as 
tabulations  of  these  parameters. 

Reduced  data  were  plotted  as  curves  of  the  parameters  versus  time, 
and  are  appropriately  grouped  in  Appendix  C  to  assist  in  comparison  of 
duplicate  data  or  of  corresponding  parameters.  Overlapping  of  curves 
plotted  on  a  single  axis  system  for  comparison  has  been  minimized  by 
use  of  offset  ordinate  scales. 

Specific  information  from  the  parameter-time  curves  has  been 
tabulated.  Data  in  the  tables,  as  well  as  in  the  curves,  are  identi¬ 
fied  by  the  symbols  used  in  Figs.  2.1  and  2.2  and  defined  in  a  footnote 
to  Table  2.6.  However,  the  tabular  information  is  arranged  to  empha¬ 
size  comparable  information  rather  than  in  the  order  position  in  the 
gage  arrays. 

2.4,1  Shot  1 


Data  from  all  19  information  channels  at  Sta  3-285  for  Shot  1 
are  presented  in  Tables  2.1  and  2.2  and  in  Figs.  C.l  through  C*7» 
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Fig.  2.4  Three -Component  Accelerometer  Mount 
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TABLE  2.1  -  Free-Field  Earth  Stress  and.  Strain  Data  -  Shot  1 
(Ground  range:  1450  ft) 


mmmmam 

Arrival 

Rise 

Peak 

Peak 

Positive  phase 

time 

time 

stress 

strain 

duration 

Gage* 

Componentt 

■sn 

Kffl 

BiH 

(sec) 

PG-9 

Air 

0.2810 

0.016 

14.45 

0.58 

pressure 

EP-1 

Vert 

0.2910 

0.019 

31.56 

_ 

0.64 

ES-1 

Vert 

0.2915 

0.026 

- 

6.0 

0.71 

EP-7 

Vert 

0.2962 

0.014 

27.43 

- 

0.62 

ES-7 

Vert 

0.2984 

0.020 

- 

17.0 

0.76 

EP-2 

Rad 

0.2931 

0.019 

9.23 

- 

O.69 

ES-2 

Rad 

0.2944 

0.034 

- 

2.7 

0.71 

EP-3 

VR  DT 

0.2943 

0.020 

17.29 

- 

O.63 

es-6 

VR  DT 

0.2990 

0.01 2 

10.0 

0.66 

EP-4 

Tang 

0.2946 

0.020 

6.96 

- 

0.62 

ES-4 

Tang 

0.2977 

0.004 

- 

1.3 

0.6l 

EP-8 

Tang 

0.2955 

0.014 

9.54 

- 

O.63 

EP-5 

TR 

0.2953 

0.007 

8.98 

- 

0.64 

ES-5 

TR 

0.2979 

0.005 

- 

2.2 

0.90 

EP-b 

VR  DA 

0.2966 

0.012 

18.30 

- 

0.66 

ES-3 

VR  DA 

0.2946 

0.011 

- 

6.6 

O.91 

TABLE  2.2  -  Pree-Fleld  Earth  Acceleration  Data  -  Shot  1 
(Ground  range :  lhJO  ft) 


■ 

Direction^ 

Arrival 

time 

Peak 

acceleration 

Peak 

velocity 

Displacement 

Peak 

Residual] 

mumu 

(g) 

KOH1 

ufsaui 

EAV 

up 

down 

0,2872 

1.51 

2.31 

m 

0.82 

O.56 

0.24 

EAR 

in 

out 

0,2872 

0.70 

1.37 

0.42 

0.26 

0.33 

0.30 

0.05 

EAT 

CCW 

CW 

0.2877 

O.63 

0.31 

0.18 

0.22 

0.04 

*  See  Fig.  2.1 
t  See  footnote  Table  2.6 
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I 


2.4.2  Shot  9 

Only  21  of  the  24  end  instruments  at  Sta  F-281  were  operational 
during  Shot  9  because  their  recorder  channels  were  assigned  to  other 
instrumentation.  Consequently  no  data  were  obtained  from  earth  Btress 
gages  EP-1,  EP-8,  and  EP-11.  Circuit  failures  external  to  the  gages 
rendered  two  of  the  remaining  channels  inoperative,  and  no  data  were 
obtained  from  stress  gage  EP-5  and  strain  gage  ES-5.  Parameter-time 
curves  for  the  19  channels  which  yielded  information  are  included  in 
Pigs.  C.8  through  C.14,  and  corresponding  numerical  data  appear  in 
Tables  2.3  and  2.4. 


TABLE  2.3  -  Free-Field  Earth  Stress  and  Strain  Data  -  Shot  9 
(Ground  Range:  1430  ft) 


Arrival  Rise  Peak 
time  time  stress 


Peak 

strain 


Positive  phase 
duration 


e*  Component 


Air 

pressure 

Air 

pressure 


EP-1 

EP-2 

ES-1 

EP-10 

ES-7 

EP-3 

ES-2 

EP-4 

ES-3 

EP-5 

ES-4 

EP-6 

ES-5 

EP-7 

es-6 

EP-8 

EP-9 

EP-11 

EP-12 


Vert 
Vert 
Vert 
Vert 
Vert 
Rad 
Rad 
TO  DT 
TO  DT 
Tang 
Tang 
TR  OCW 
TR  OCW 
TO  DA 
TO  DA 
TO  DA 
TO  OCCW 
TV  DCCW 
TV  DCW 


Not  connected  to  recorder 


1.2967 

0.012 

19.46 

1.3011 

0.018 

- 

1.2949 

0.013 

19.49 

1.2952 

0.012 

- 

1.3000 

0.005  ; 

3.60 

1.3022 

0.003 

- 

1.2911 

o.oi4  ; 

8.64 

1.2911 

0.013 

- 

Circuit  failure  I 

0.002 

1.2983  I  0.005  I  3.42 
Circuit  failure 
1.2964  |  0.006  I  11.86 
1.2975  I  0.005  I  - 

Not  connected  to  recorder 
1.2959  !  0.005  !  4.64  | 

Not  connected  to  recorder 
0.007 


*  See  Fig.  2.2 
t  See  footnote  Table  2.6 
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TABLE  2=4  -  Free-Fiel d  Earth  Acceleration  Bata  -  Shot  9 
(Ground  Range:  1430  ft) 


Gage* 

— — 

Directiont 

Arrival 

time 

Peak 

acceleration 

Peak 

velocity 

Peak 

Residual  I 

uxmu 

mow 

16391 

■earn 

EAV 

up 

down 

1.2935 

0.70 

4.34 

o.n 

0.55 

0.00 

0.57 

0.25 

EAR 

in 

out 

1.2933 

0.80 

1.65 

33 
•  • 

O  O 

0.22 

0.03 

0.02 

EAT 

CCW 

CW 

1.2935 

0.57 

0.46 

0.08 

0.06 

0.08 

0.05 

0.05 

*  See  Fig.  2.2 
t  See  footnote  Table  2.6 


2.4.3  Shot  10 


All  24  information  channels  at  Sta  F-281  were  operating  during 
Shot  10  and  produced  valid  data,  given  in  Figs.  C.15  through  C.22  and 
Tables  2.5  and  2.6. 

Partition  of  the  air  overpressure  from  Shot  10  into  a  precursor 
shock  and  main  shock  as  indicated  in  Fig.  C.15  required  subdivision  of 
arrival  times  and  peak  values  of  the  parameters  (Tables  2.5  and  2.6). 
The  positive  phase  durations  in  Table  2.5  represent  the  time  interval 
from  arrival  of  the  precursor  to  crossover  into  the  negative  phase  of 
the  phenomenon. 

2 . 5  DISCUSSION  OF  RESULTS 


Analytical  operations  to  which  measured  parameters  must  be  sub¬ 
jected  to  evaluate  their  usefulness  are  defined  in  general  terms  in 
Section  2.2.  Comparison  of  incident  overpressure  and  vertical  stress 
measurements  will  verify  the  source  of  observed  earth  parameters.  A 
measure  of  reliability  of  the  instrumentation  as  placed  can  be  obtained 
from  comparison  of  data  from  duplicate  component  gages.  Comparisons  of 
corresponding  earth  stress  and  strain  components  can  indicate  the  elasto- 
plastic  behavior  of  the  soil. 

Derivation  of  principal  stresses  was  feasible  for  only  one  of  the 
three  sets  of  data  because  insufficient  component  measurements  were  ob¬ 
tained  from  Shots  1  and  9«  The  principal  stress  can  be  approximated 
for  Shot  1  if  radial  symmetry  is  assumed,  but  the  angle  of  34°  between 
the  blast  line,  which  controlled  orientation  of  the  earth  gages  for 
Shot  9»  end  the  line  from  true  ground,  zero  through  the  center  of  the 
array  precludes  computation  of  the  principal  stresses  (Fig.  1.2). 
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See  Fig.  2.2. 

See  footnote  Table  2.6 


TABLE  2.6  -  Free-Field  Earth  Acceleration  Data  -  Shot  10  (Ground  range:  1120  ft) 
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Comparison  of  air  overpressure  and  vertical  earth  stress  for  each 
shot  shows  that  these  parameters  vary  in  essentially  similar  patterns, 
although  the  incident  overpressures  are  lower.  There  is  also  close 
similarity  in  the  time  variations  of  most  of  the  other  stress  and  strain 
components,  as  evidenced  by  the  curves  in  Appendix  C.  These  facts  es¬ 
tablish  conclusively  that  all  but  minor  portions  of  the  earth  stresses 
and  strains  resulted  directly  from  incidence  of  the  air  shock  wave  at 
the  ground  surface  adjacent  to  the  gage  array. 


2.5.1  Reliability  of  Duplicate  Observations 

Comparison  of  duplicate  observations  of  the  same  parameters  as 
a  function  of  time  is  one  measure  of  the  ability  of  similar  gages, 
placed  under  similar  conditions  and  subject  to  similar  loading  histories, 
to  register  similar  parame ter-time  patterns.  Variations  of  individual 
gage  characteristics  are  ignored,  since  reliability  of  the  gages  under 
repeated  static  calibration  and  transient  unloading  tests  is  probably 
better  than  97  per  cent. 

The  effect  of  the  presence  of  gages  and  of  disturbance  of  a 
volume  of  soil  inherent  to  installing  them  is  also  ignored.  It  is  real¬ 
ized  that  under  certain  conditions  these  factors  may  so  distort  the 
stress  distribution  in  the  soil  as  to  cause  the  gages  to  indicate 
stresses  differing  by  factors  as  high  as  2  or  3  from  those  in  the  un¬ 
disturbed  medium.  However,  magnitudes  of  these  effects  have  been  meas¬ 
ured  and  computed  for  static  loads  only,  and  projection  of  the  assump¬ 
tions  on  which  such  estimates  are  based  to  distortion  of  transient 
stress  distribution  may  be  questionable  in  the  present  state  of  knowl¬ 
edge  of  dynamic  behavior  of  soils. 

The  instrumentation  plan  for  Shot  1  (Fig.  2.1)  included  three 
sets  of  duplicate  observations:  EP-1  and  BP-7,  spaced  18  ft  apart, 
measured  vertical  earth  stress;  EP-4  and  EP-8,  12  ft  apart,  measured 
tangential  stress;  and  ES-1  and  ES-7  measured  vertical  earth  strain. 
Duplicate  observations  were  limited  for  Shot  9  (Fig.  2.2)  to  vertical 
stress,  EP-2  and  EP-10,  and  vertical  strain,  ES-1  and  ES-7.  Data  from 
Shot  10  (Fig.  2.2)  provided  more  comparable  observations :  three  ver¬ 
tical  stresses,  EP-1,  EP-2,  and  EP-10;  two  45°  vertical-radial  stresses, 
EP-7  and  EP-8;  and  two  vertical  strains,  ES--1  and  ES-7. 

The  study  of  reliability  is  based  on  comparison  of  data  samples 
from  duplicate  gages.  Sample  data  were  taken  throughout  the  positive 
phase  of  measured  stresses  and  strains  and  were  chosen  by  identity  of 
shape  of  minor  features  of  the  record  trace  pairs  as  well  as  by  approxi¬ 
mate  simultaneity.  Mean  values  of  sample  pairs  and  deviations  of  the 
data  from  these  means  were  computed.  Reproducibility  of  the  data  is 
expressed  as  the  average  of  the  deviations  and  as  the  range  of  devia¬ 
tions.  Results  of  this  study  are  presented  in  Table  2.7.  The  only 
features  of  the  data  in  Table  2.7  which  require  clarification  are  the 
effects  of  residual  strains,  and  the  possibly  anomalous  stresses  indi¬ 
cated  by  EP-1  for  Shot  10. 
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TABLE  2.7  -  Comparison  of  Duplicate  Data 


.  1 

Parameter 

Number  of 

Deviation  from  mean  (%) 

Shot 

component 

observations 

Range 

Average 

1 

Vert  stress 

72 

4-42 

16 

Tang  stress 

67 

o-31 

13 

Vert  strain 

66 

10-54 

35 

9 

Vert  stress 

20 

10-22 

6 

Vert  strain 

15 

16-32 

26 

10 

Vert  stress* 

13 

25-91 

55 

Vert  stress t 

13 

0-18 

8 

Vert-rad  stress 

14 

1-36 

14 

Vert  strain 

14 

3-21 

14 

*  Data  from  EP-1  and  EP-2 
t  Data  from  EP-2  and  EP-10  (see  text) 


Extension  of  the  strain-time  curves  beyond  the  negative  phase 
shows,  in  all  cases,  a  residual  strain,  probably  indicative  of  plastic 
defoliation  of  the  soil  in  the  vicinity  of  the  gages.  It  is  also  noted 
in  Appendix  B  that  placement  of  earth  strain  gages  to  respond  to  the 
vortical  component  involves  considerable  difficulty  in  attaining  uni¬ 
formity  of  compaction  of  the  soil  between  the  end  plates  while  retain¬ 
ing  tile  proper  mechanical  setting  of  the  transducer.  Consequently 
performance  of  two  such  gages  in  response  to  approximately  equal  tran¬ 
sient  loads  may  indicate  a  much  wider  spread  of  strains  than  would  be 
anticipated  from  corresponding  stress  measurements.  Adjustment  of  the 
vertical  strain  data  from  duplicate  gages  to  equal  residual  strains 
makes  the  deviation  of  the  data  from  the  mean  much  less  than  it  is  for 
the  raw  data.  ES-1  and  ES-7  showed  residual  strains  for  Shot  1  of 
0.4  x  10  '  and  5  x  10"^  in. /in.,  and  the  peak  strain  after  correction 
was  twice  as  great  for  ES-7  as  for  ES-1.  These  data,  as  well  as  the 
consistent  one-sidedness  (ES-7  always  greater)  of  the  entire  strain¬ 
time  comparison,  suggest  a  generally  lower  compacted  density  of  the 
soil  in  the  vicinity  of  ES-7  and  particularly  beneath  the  upper  end 
plate  of  that  gage.  Residual  strains  from  the  vertical  gages  for  Shots 

9  and  10  were  more  comparable:  they  were  about  5  x  10~4  in. /in.  for 
Shot  9  with  a  difference  of  only  10”^  in. /in.  between  ES-1  and  ES-7 
and  8  x  10”^  in. /in.  for  Shot  10  with  approximately  the  same  difference 
between  readings. 

Two  of  the  vertical  stress  gages,  EP-1  and  EPr2,  at  Sta  F-281 
(Fig.  2.2)  were  separated  by  only  3  ft  while  the  third,  EP-10,  was  24 
ft  from  EP-2.  However,  EP-2  and  EP-10  each  showed  stresses  for  Shot 

10  which  were  about  1.7  times  those  indicated  by  EP-1  acid  1.5  times 


42 

SECRET  -  RESTRICTED  DATA 


the  incident  air  overpressure .  The  apparent  anomaly  of  the  EP-1  stress 
data  makes  evaluation  of  information  from  the  three  gages  difficult 
since  labeling  the  EP-1  data  anomalous  may  be  questioned  in  the  light 
of  close  agreement  vith  incident  air  overpressure.  However,  vertical 
earth  stresses  at  5  ft  exceeded  incident  air  overpressures  by  20  to  100 
per  cent  for  Shots  1,  9,  and  10  of  this  series  and  Shot  7  of  TUMBLER- 
SNAPPER.  There  is  some  reason  to  expect  that  vertical  earth  stress 
might  exceed  incident  overpressures  by  as  much  as  100  per  cent  because 
of  differences  in  acoustic  impedance  across  the  air-earth  boundary. 
Principal  stresses  computed  from  Shot  10  data  have  also  given  strong 
support  to  the  higher  values  of  vertical  earth  stress.  Consequently 
data  from  EP-2  and  EP-10  are  considered  the  most  reliable. 

2.5.2  Stress -Strain  Relations 


Curves  of  corresponding  components  of  earth  stress  as  a  function 
of  strain  give  information  concerning  the  elastic  and  plastic  response 
of  soils  to  loading.  Plastic  deformation  of  the  soil  under  load  im¬ 
plies  time  lags  between  application  of  a  stress  and  stabilization  of 
the  corresponding  strain.  This  lag  is  taken  into  account  in  measuring 
static  stress-strain  characteristics  by  retaining  a  series  of  incre¬ 
mental  loads  for  periods  that  vary  from  minutes  to  days  and  result  in 
approximately  complete  equilibrium  for  each  increment.  Transient  loads 
produced  in  soil  by  shock  obviously  do  not  permit  stabilization  of 
plastic  strain  because  they  change  too  rapidly.  Consequently,  meas¬ 
ured  transient  earth  strains  may  be  expected  to  lag  corresponding  meas¬ 
ured  stresses.  These  lags  are  random  and  small  and  usually  not  much 
greater  tlian  differences  in  arrival  times.  Average  measured  deviations 
of  the  strain  lags  from  arrival  time  differences  were  about  1.3  msec. 

It  is  practicable  to  correlate  major  and  minor  features  of  correspond¬ 
ing  stress-  and  strain-time  curves  such  as  peaks  and  some  features  of 
the  precursor  shocks,  disregarding  exact  time  correlations.  This  type 
of  correlation  is  not  dependable  within  the  sharply  rising  front  of 
the  shock  wave  or  in  the  smoother  portion  of  a  clean  wave  beyond  the 
peak. 

Uncertainty  of  correlation  makes  doubtful  the  true  shape  of 
the  loading  branch  (initial  loading  to  peak)  of  a  stress-strain  curve 
as  well  as  some  portions  of  the  unloading  branch  (peak  to  negative 
phase).  Consequently  such  regions  of  the  curves  are  represented  by 
straight  lines.  The  curvature  and  minor  loops  of  the  precursor  shock 
region  (Shot  10)  and  the  unloading  branch  represent  spikes  and  troughs 
and  are  reasonably  dependable.  Fortunately  most  of  the  stress-  and 
strain-time  curves  were  sufficiently  ragged  beyond  the  peaks  and 
throughout  the  precursor  shock  to  facilitate  correlation.  The  stress- 
strain  curves  form  open  hysteresis  loops  in  which  the  enclosed  areas 
between  loading  and  unloading  branches  represent  energy  dissipated 
through  deformation  of  the  soil.  The  loops  are  not  closed  because 
plastic  deformation  resulted  in  residual  strains.  Arbitrary  closure 
of  the  loops  permits  estimation  of  the  energy  dissipated  through  hys¬ 
teresis. 

Stress-strain  curves  were  plotted  for  a  single  set  of  vertical 
component  measurements  for  Shot  1  (Fig.  2.5)  and  for  duplicate  vertical 
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Fig.  2.5 


Vertical  Earth  Stress 
vs  Strain  -  Shot  1 


and  orthogonal  vertical-radial  obli¬ 
que  data  for  both  Shots  9  and  10. 

All  of  these  curves  show  hysteresis 
and  residual  strain  (circled  end 
points  of  solid  curves).  Equal  areas 
within  the  several  hysteresis  loops 
as  drawn  do  not  represent  similar 
stress-strain  energies  because  of 
scaled  differences.  The  product  of 
the  units  in  which  the  stresses  and 
strains  are  presented  (psi  and  10“^ 
in. /in.)  result  in  energy  density 
rather  than  energy,  and  areas  within 
the  hysteresis  loops  represent  micro- 
foot-pounds  per  cubic  foot. 

The  forms  of  the  stress-strain 
curves  from  vertical  components  for 
Shot  9  (Fig.  2.6)  are  similar  as  are 
those  from  the  two  sets  of  vertical- 
radial  oblique  data  (Figs.  2.7  and 
2.8).  These  curves  are  considerably 
smoother  than  that  from  Shot  1  (Fig. 
2.5)  or  those  from  Shot  10  because 
the  incident  shock  and  induced  earth 
stresses  were  relatively  clean  waves. 

Evidence  of  the  precursor 
shock  appears  in  the  stress-strain  curves  from  Shot  10  as  a  lobe  and 
some  ragged  oscillations  in  the  loading  branches.  The  duplicate  vertical 
component  curves  (Fig.  2.9)  and  the  orthogonal  oblique  pair  (Figs.  2.10 
and  2.1l)  all  exhibit  generally  similar  shapes  as  they  did  for  Shot  9. 

Energy  dissipation  densities  derived  from  the  stress-strain  curves 
are  presented  in  Table  2.8  with  corresponding  values  of  peak  stress, 
strain,  and  positive-phase  durations.  Duplicate  vertical  component  data 
arc  identified  by  the  pertinent  stress  gage  symbol.  The  energy  densities 
computed  from  the  hysteresis  curve  areas  are  probably  subject  to  errors 
of  the  order  of  50  per  cent.  The  principal  sources  of  these  errors  are 
the  uncertainties  attached  to  gage  placement  and  the  assumption  of  linear 
stress-strain  relationships  between  observed  points  on  the  curves.  The 
latter  assumption  is  especially  questionable  within  the  loading  branch 
of  the  curve  where  consolidation  or  shearing  of  the  soil  during  loading 
would  require  nonlinearity  between  stress  and  strain.  Consideration  of 
the  stress-strain  ratios  at  maximum  stress  and  the  computed  energy  den¬ 
sities  suggests  that  for  the  duplicate  vertical  component  data  the  error 
in  curve  areas  may  be  more  serious  than  in  the  measured  parameters. 
Comparison  of  the  stress-strain  ratios  for  Shot  9  shows  that  that  from 
EP-2  data  is  about  1.5  times  that  from  EP-10,  and  the  same  comparison 
for  Shot  10  gives  the  same  factor.  On  the  other  hand,  comparison  of  the 
corresponding  energy  densities  for  Shot  9  indicates  aratio  of  0.58  and 
a  ratio  of  Q.6l  for  similar  comparison  of  Shot  10  energy  densities, 
roughly  the  inverse  of  the  stress-strain  ratios. 

However,  the  energy  data  give  some  reasonable  indication  of  the 
expected  dissipation  caused  by  plastic  deformation  of  the  soil. 


44 

SECRET- RESTRICTED  DATA 


Knii  (1C r4  ia./lM.) 


Strain  (KT*  m./ln.) 

Fig.  2.6  Duplicate  Vertical  Earth  Stress  vs  Strain  -  Shot  9 


|»*»»)  9tr.ii  Utr*  IM./lM.l 


Fig.  2.7  Oblique  Earth  Stress  Fig.  2.8  Oblique  Earth  Stress 

vs  Strain  -  Shot  9  •  vs  Strain  -  Shot  9 
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TABLE  2.8  -  Stress -Strain  Energy  Dissipation 


T 


Shot 

Component 

Peak 

stress 

Peak 

strain 

Pos  phene 
duration 

Density  of 
energy  dissipa¬ 
tion 

R1SI 

nsesran 

■SH! 

IBSSBlBil 

■ 

Vert 

31.6 

6.0 

0.67 

12 

s 

Vert  ( EP-2 ) 

19.5 

7.2 

0.70 

32 

Vert  ( EP-10 ) 

19.5 

11.0 

0.68 

55 

■ 

Vert -Rad  (DT) 

8.6 

2.7 

0.92 

6.4 

■ 

Vert-Rad  (DA) 

11.9 

8.9 

O.50 

28 

10 

Vert  (EP-2) 

56.1 

13.3 

0.43 

200 

Vert  (EP-10) 

59-9 

21.2 

0.43 

330 

Vert-Rad  (DT) 

10.7 

5.2 

0.42 

l6 

Vert-Rad  (DA) 

20.1 

6.1 

0.31 

38 

Stress-  and  strain-time  curves  for  the  components  used  in  com¬ 
paring  stress-strain  relations  (Figs.  C.2,  C.9,  C.10,  C.l6,  and  C.l8) 
show  that  in  all  cases  the  transient  loads  resulted  in  permanent  com¬ 
pressive  strains  (circled  points  on  the  curves)  which  are  interpreted 
as  plastic  deformations.  These  appear  to  he  a  relatively  large  part  of 
the  peak  strains.  However,  both  peak  and  residual  strains  are  small 
when  the  magnitudes  of  the  corresponding  applied  loads  are  considered. 
Evidently  the  peak  strains  are  relatively  small  because  of  lag  in  re¬ 
sponse  of  the  soil  to  sudden  loading.  They  probably  represent  mostly 
elastic  strain,  which  is  recovered  before  the  end  of  the  positive  phase 
of  the  transient  stress  but  to  which  is  added  throughout  the  positive 
phase  plastic  strain  which  is  not  recovered  except  as  the  low-intensity 
forces  of  the  negative  phase  may  affect  it.  Lag  in  deformation  of  the 
soil  under  the  transient  load  accounts  for  the  appreciably  longer  posi¬ 
tive-phase  durations  for  strain  gage  data. 

2.5-3  Stress  Tensors 


Stress  measurements  from  Shot  1  included  components  for  deter¬ 
mining  the  principal  stresses  if  radial  symmetry  is  assumed.  The  meas¬ 
ured  tangential  earth  stress  is  considered  to  be  the  minor  principal 
stress,  and  data  from  the  vertical,  radial,  and  vertical-radial  oblique 
stress  gages  are  used  as  though  a  plane  stress  system  existed,  ,  Compu¬ 
tations  are  similar  to  those  described  in  an  earlier  report. Re¬ 
sults  of  this  analysis  showed  that  the  major  principal  stress  was  only 
slightly  greater  than  the  mean  measured  vertical  stress,  and  that  it 
was  directed  between  10°  and  20°  from  the  vertical.  These  results  are 
based  on  data  read  from  the  stress-time  curves  at  identical  times  and 
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asymmetry  of  the  shock  front. 

Identity  of  the  stress -time  pattern,  as  indicated  by  similarity 
of  peak  or  trough  shape  and  approximate  position  in  the  stress -time 
sequence,  defined  comparable  stresses  for  different  components  of  the 
Shot  10  data.  Data  representing  a  single  nominal  instant  actually  in¬ 
volved  a  range  of  times  as  great  as  10  msec  over  the  group  of  stress - 
time  curves.  However,  in  view  of  the  spread  of  the  gage  array  and  the 
evident  correspondence  of  stress -time  patterns  for  different  gages,  the 
assumption  of  simultaneity  was  considered  reasonable,  and  mean  values 
were  used  to  define  the  nominal  time  for  each  set  of  stress  data.  Only 
events  which  could  be  readily  identified  on  all  or  nearly  all  stress¬ 
time  curves  (Figs.  C.15  through  C.20)  were  considered  suitable  for  anal¬ 
ysis.  Data  representing  14-  such  events  during  the  positive  phase  of 
stress -time  curves  for  Shot  10  were  compiled  for  use  in  computing  prin¬ 
cipal  stresses. 

The  theoretical  background  and  procedure  for  deriving  principal 
stresses  for  three-dimensional  conditions  are  described  in  Appendix  D; 
definitions  of  the  stress  system  are  included  there  to  refresh  the 
reader's  memory  of  basic  stress  relations.  The  basis  for  these  rela¬ 
tions  is  the  existence  of  stress  equilibrium.  Equilibrium  i&  not  nec¬ 
essarily  a  general  condition  during  transient  loading  of  a  soil  mass, 
but  it  probably  does  exist  instantaneously  throughout  the  load  cycle 
and  is  probably  more  fully  realized  at  maxima  and  minima  than  during 
intermediate  portions  of  the  cycle.  Therefore  the  assumption  that  the 
derived  data  represent  conditions  when  all  stresses  are  in  equilibrium 
is  reasonable  and  permits  analysis  as  described  in  Appendix  D. 

Data  for  stress  analysis  should  satisfy  tests  for  the  static 
equilibrium  of  stresses.  The  simplest  test  for  such  equilibrium  is 
equality  of  the  suras  of  the  normal  stresses  on  all  triads  of  mutually 
orthogonal  planes  through  a  point.  This  is  only  one  index  and  is  a 
necessary  but  not  a  sufficient  condition  for  equilibrium.  Its  sim¬ 
plicity  and  the  fact  that  it  may  be  expected  to  yield  indicative,  If 
not  conclusive,  evidence  of  equilibrium  recommended  its  use. 

Observed  stresses  must  accordingly  satisfy  the  equality 

aV  +  °R  +  ^T  =  aVR  DA  +  °VR  DT  +  aT 

=  aVT  DCW  +  aVT  DCCW  +  aR 
=  aRT  OCW  +  aRT  OCCW  +  V 

or  the  more  practical  forms  in  which  common  terms  in  each  pair  of  sums 
are  eliminated: 

(ffV  +  °r)  "  (°VR  DA  +  °W  Dt)  =  A0VR  =  0 
(°V  +  CTT  )  "  (aVT  DCW  +  °VT  DCCw)  =  ^VT  3  0 
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(CR  +  ctt)  "  (°RT  OCW  +  aRT  OGCwj  =  *°RT  “  °* 

It  ia  recognized  that  anticipated  instrumentation  errors  might  he  domi¬ 
nant  and  yield  inequalities.  Random  inequalities  of  magnitudes  similar 
to  those  of  anticipated  errors  could  he  regarded  as  indicative  of  ap¬ 
proximate  stress  equilibrium.  However,  absence  of  randomness  must  be 
regarded  as  principally  indicative  of  nonunifom  stress  distribution 
and  can  not  distinguish  between  effects  of  gage  placement  and  nonequi- 
librium  of  stresses. 

Inequalities  AtJyg,  Ac^,,  and  La^j>  are  derived  for  each  of  the 
14  sets  of  observations .  These  differences  are  normalized  to  give  com¬ 
parable  values  of  La  f  Sc  where  So  represents  the  arithmetic  mean  of  the 
summed  stress  pairs,  i.e., 

pV  +  <*R  +  °VR  DA  +  ffVR  Dt| 

2ara  =[ - 2 - y 

which  may  be  thought  of  as  the  sum  of  each  stress  pair  if  streBB  equi¬ 
librium  were  realized.  The  mean  of  these  normalized  differences  and 
the  probable  deviation,  e,  of  any  observation  from  the  mean  are  the 
criteria  of  equilibrium.  Three  sets  of  vertical  stress  data  -  EP-1 
alone,  the  mean  of  EP-1,  EP-2,  and  EP-10,  and  the  mean  of  EP-2  and  EP-10 
data  (of  Fig.  2.2  and  discussion  of  Shot  10  data  in  Section  2.5*1)  - 
are  used  for  the  equilibrium  tests. 

Results  of  this  analysis  are  tabulated  (Table  2,9)  to  compare 
mean  values  of  normalized  inequalities  with  probable  deviations  of  any 
of  these  observations  from  the  mean.  Differences  in  BtresB  sums  for 
EP-1  and  (TRm  data  are  fairly  random,  those  for  the  latter  being  nearly 
evenly  balanced  between  positive  and  negative  differences.  Differences 
derived  from  averaged  vertical  stress  observations  are,  on  the  contrary, 
very  one-sided,  sums  including  cjy  being  greater  in  all  instances. 


TABLE  2.9  -  Stress  Equilibrium  Check  -  Shot  10 


Source  of  data 

I'll 

eVR 

AgVT 

eVT 

eRT 

HKJJH 

■OB 

BTC9JB 

StlJBB 

10] 

EP-1 

(EP-1  +  EP-2  +  EP-10 )/2 
(EP-2  +  EP-10 )/2 

34 

4  6 

68 

21 

19 

19 

30 

47 

70 

32 

11 

7 

■  ■  .  , 

6 

6 

6 

18 

18 

18 
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Table  2.9  shows  no  conclusive  evidence  either  for  or  against 
the  stress  equilibrium.  On  the  basis  of  the  mean  of  the  normalized 
difference,  which  is  sensitive  to  randomness  of  sign,  data  from  EP-1 
show  departures  from  equilibrium  of  the  order  of  the  anticipated  instru¬ 
mentation  error.  Probable  deviations  from  the  mean  normalized  differ¬ 
ences  are  insensitive  to  sign,  and  the  third  set  of  data  is  considerably 
better  than  those  for  EP-1  and  somewhat  better  than  mean  data  from  all 
three  gages.  But  although  these  values  indicate  deviations  well  within 
the  range  of  anticipated  instrumentation  error,  they  can  not  distinguish 
the  equilibrium  criterion  from  either  instrumentation  error  or  distorted 
stress  distribution  in  the  vicinity  of  instrumentation.  Perturbation 
of  stress  distribution  caused  by  recompaction  of  the  soil  mass  during 
gage  installation  or  the  presence  of  gages  within  the  soil  may  result 
in  relatively  large  one-sided  deviations  of  stress  sums  without  inter¬ 
fering  with  instantaneous  stress  equilibrium. 

Principal  stresses  representative  of  the  stress  tensor  in  the 
earth  at  14  instants  after  arrival  of  the  air  shock  above  the  gage  ar¬ 
ray  were  computed  despite  the  indeterminate  results  of  the  stress  equi¬ 
librium  check.  These  computations,  outlined  in  Appendix  D,  include 
derivation  of  (l)  shear  stresses  in  the  V-,  R-,  and  T-planeB,  (2)  major, 
intermediate,  and  minor  principal  stresses,  (3)  angles  between  each 
principal  stress  and  the  three  coordinate  axes,  and  (4)  maximum  shear 
stress.  The  first  two  were  performed  for  three  sets  of  vertical  stress 
data:  EP-1  data  only;  the  mean  of  EP-1,  EP-2,  and  EP-10  data;  and  the 
mean  of  EP-2  and  EP-10  data.  The  remaining  computations  were  made  using 
only  the  mean  data  from  EP-2  and  EP-10. 

Shear  stresses  on  the  vertical,  radial,  and  tangential  planes, 
derived  from  measured  normal  stresses  as  described  in  Appendix  D,  are 
shown  as  functions  of  time  in  Pig.  2.12.  Stress  equilibrium  (assumed 
if  not  rigorously  proven)  requires  the  equalities  typ  =  tKy,  Tyj,  = 
and  =  ir^.  These  shear  stresses  and  the  normal  stresses,  ay,  aR, 
and  c*j>  appear  in  the  cubic  equation  (Eq  D.5)  from  which  the  principal 
stresses  are  derived  (Appendix  D). 

Analysis  based  on  cjy  data  from  EP-1  resulted  in  imaginary  roots 
of  the  cubic  equation  at  one  or  more  points  near  the  end  of  the  positive 
phase  and  in  negative  stresses  (tensions)  during  the  precursor  portion 
of  the  minor  principal  stress  as  great  as  85  per  cent  of  the  maximum 
compressive  stress  in  that  part  of  the  shock  wave.  Shear  stresses  in 
the  vertical  plane  were  large,  having  an  average  value  nearly  equal  to 
cry.  Derived  angles  between  the  major  principal  stress  and  the  axes 
were  erratic,  and  except  in  the  Immediate  vicinity  of  the  main  shock 
peak,  shoved  the  stress  to  be  inclined  almost  6 0°  to  the  vertical  axis. 
The  major  principal  stress  differed  radically  from  cy  in  magnitude  as 
well  as  direction.  Consequently  if  the  reasoning  which  recommended 
performance  of  this  experiment  in  connection  with  an  air  burst  weapon 
is  sound,  either  EP-1  data  lead  to  false  results  -  a  likely  conclusion 
in  view  of  the  tensions,  imaginary  roots,  high  shear  stresses  in  the 
V-plane,  and  departure  of  the  major  principal  stress  from  vertical  - 
or  induced  stresses  are  propagated  through  the  earth  near  the  surface 
at  a  relatively  steep  angle  despite  generation  at  the  surface. 
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Choice  between  data  representing  means  for  the  three  vertical 
gages  end  those  derived,  from  EP-2  and  iSP-10  alone  is  not  so  obvious  as 
was  elimination  of  EP-1  data.  Average  shear  stresses  in  the  vertical 
plane  were  somewhat  smaller  for  the  2-gage  mean  data,  and  probable  de¬ 
viation  of  the  normalized  differences  of  stress  sums  in  Table  2.6  was 
also  more  favorable  for  them.  Directions  of  the  major  principal  stresses 
derived  from  the  two  sets  of  data  differed  by  only  about  5°  although 
the  3-gage  data  gave  the  larger  inclination  to  the  vertical  axis.  There¬ 
fore  a  conclusive  choice  is  not  possible,  but  the  2-gage  mean  data  seem 
the  better. 

Curves  for  principal  stresses  derived  from  mean  data  from  EP-2 
and  EP-10  as  a  function  of  time  (Figs.  2.13,  2.14,  and  2.15)  are  self- 
explanatory.  Similar  curves  for  the  angle  between  stresses  and  coordi¬ 
nate  axes  in  the  same  figures  need  discussion  only  from  the  standpoint 
that  orientation  of  the  two  lesser  principal  stresses  is  less  precise. 
Deriving  the  angles  in  the  latter  cases  involves  operations  with  large 
differences  of  small  numbers,  and  results  are  much  more  sensitive  to 
inaccuracies  in  derived  stresses.  However,  these  stresses  and  their 
directions  are  relatively  insignificant  because  the  major  principal 
stress  is  greater  than  either  by  a  factor  of  at  least  0.  The  inter¬ 
mediate  and  minor  principal  stresses  are  thus  essentially  of  only  aca¬ 
demic  interest  under  the  conditions  produced  by  air  shock  loading  on 
the  earth. 

It  is  of  interest  to  note  the  differences  in  general  direction 
of  the  two  lesser  principal  stresses  during  the  precursor  and  main 
shock  phases  of  loading.  During  the  precursor  shock  the  intermediate 
stress  is  directed  nearly  midway  between  the  radial  and  tangential  axes; 
during  the  main  shock  it  is  very  nearly  parallel  to  the  tangential  axis. 
The  minor  principal  stress  shows  a  similar  trend:  nearly  parallel  to 
the  radial  axis  during  the  main  shock,  between  the  radial  end  tangen¬ 
tial  axes  during  the  precursor  shock.  Inaccuracies  of  computations  for 
precursor  data  are  sufficient  to  account  for  the  erratic  oscillation  of 
the  angles  for  cr2  and  but  it  seems  pi-obable  that  trends  described 
are  real  if  not  strongly  defined. 

Shear  stresses  derived  from  principal  stresses  in  accordance 
with  the  equation 

_  K  -  c3} 

^raax  2 

are  maxima  and  act  on  the  plane  which  bisects  the  angle  between  major 
and  minor  principal  planes.  The  maximum  shear  stress,  i,nQX,  is  shown 
as  a  function  of  time  in  Fig.  2.16. 

Magnitudes  of  the  maximum  shear  stress  are  several  times  great¬ 
er  than  those  of  the  shears  derived  from  measured  normal  stresses  (Fig 
2.12).  This  situation  is  contrary  to  that  between  the  major  principal 
stress  and  the  measured  normal  stress  on  the  V-plane.  However,  consi¬ 
deration  of  the  directions  of  the  planes  in  which  the  derived  shears 
act  show  them  to  be  more  nearly  coincident  with  the  principal  planes  - 
in  which  shears  vanish  -  than  is  the  plane  of  maximum  shear. 
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13  Major  Principal  Stress  and  Orientation  vs  Time  -  Shot  10 


52 

SECRET- RESTRICTED  DATA 


Time  Iwrt 


Fig.  2.16  Maximum  Shear  Stress  vs  Time  -  Shot  10 


2.5.h  Ground.  Motion 


Particle  velocity  and  displacement  may  be  derived  from  measured 
particle  acceleration  by  integration  with  respect  to  time .  A  velocity¬ 
time  curve  resulting  from  integration  of  acceleration-time  data  is  sub¬ 
ject  to  the  restriction  that  the  velocity  must  ultimately  become  zero. 
Normal  field  data  frequently  require  some  adjustment  to  satisfy  the 
final  zero  velocity  requirement. 

Displacement-time  curves  are  derived  from  the  corrected  velocity¬ 
time  data  by  a  similar  integration  process.  Adjustment  of  displacements 
is  less  clear-cut  than  for  velocities.  Frequently  no  adjustment  is  nec¬ 
essary  to  provide  a  displacement-time  curve  which  appears  reasonable. 
However,  it  is  not  always  feasible  to  make  the  proper  velocity  correc¬ 
tion  from  inspection  of  that  curve  only,  and  consideration  of  first 
approximation  curves  of  both  velocity  and  displacement  may  be  necessary 
for  logical  correction.  The  simplest  correction  is  an  assumed  zero 
shift  of  the  acceleration  curve  which  will  have  the  effect  of  linear 
increase  or  decrease  in  velocity,  and  results  in  a  parabolic  correction 
of  the  displacement  curve.  Correction  may  also  be  made  by  assuming  a 
linear  zero  drift  in  the  acceleration  data  which  results  in  a  parabolic 
correction  in  the  velocity  curve  and  a  third-power  correction  in  the 
displacement  curve.  Combinations  of  these  corrections  may  be  helpful, 
but  usually  the  less  complex  corrections  yield  reasonable  results  pro¬ 
vided  the  original  data  are  reliable  and  the  time  interval  is  short. 

Integration  and  adjustment  of  the  data  may  be  ,UBed  with  reason¬ 
able  success  for  data  from  phenomena  of  a  few  seconds  duration.  How¬ 
ever,  durations  of  the  order  of  tens  of  seconds  are  difficult  to  handle 
because  errors  are  cumulative  and  real  basic  information  may  become 
hopelessly  lost  in  the  effects  of  minor  correction  errors. 
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Free -field  measurements  of  vertical,  radial,  and  tangential 
components  of  earth  acceleration  were  made  during  Shots  1,  9,  and  10. 
Parameter-time  curves  of  acceleration,  velocity,  and  displacement  for 
Shot  1  are  presented  in  Figs.  C.5  -  C.7;  for  Shot  9  in  Figs.  C.12  - 
C.l4j  and  for  Shot  10  in  Figs.  C.21  -  C.23.  Arrows  on  parameter  axes 
of  these  curves  indicate  direction  of  motion. 

Maximum  accelerations,  velocities,  and  displacements  and  resid¬ 
ual  displacements  for  each  component  are  compiled  in  Tables  2.2,  2.4, 
and  2.6.  Residual  displacements  sure  those  at  which  each  displacement 
is  assumed  to  remain  constant  since  they  occur  at  times  after  acceler¬ 
ations  and  velocities  are  assumed  to  become  vanishingly  small. 

Acceleration-time  curves  for  Shots  1  and  9  show  no  appreciable 
departures  from  zero  following  crossover  after  the  positive  phase  of 
the  stress  data.  However,  sufficient  deviation  from  zero  appears  in 
several  of  the  velocity-time  curves  beyond  the  positive  stress  phase 
to  suggest  integration  to  a  later  time. 

The  period  during  which  the  velocities  are  assumed  to  be  sig¬ 
nificant  are  represented  by  the  solid  curves  for  both  velocity  and 
displacement;  the  dashed  portions  indicate  continuing  trends  resulting 
from  long-period  minor  departures  of  the  acceleration  curves  from  zero. 
These  departures  emphasize  one  of  the  difficulties  of  properly  adjust¬ 
ing  the  acceleration  zero  to  give  reasonable  velocity  data. 

Integration  of  acceleration  data  from  Shot  10  measurements  has 
been  carried  well  beyond  the  start  of  the  negative  phase  because  a 
relatively  strong  secondary  signal  was  recorded  at  about  0.9  sec. 

This  signal  corresponds  to  the  strong  secondary  events  noted  at  about 
the  same  time  on  all  stress  and  strain  measurements,  probably  a  second¬ 
ary  shock  in  the  air  wave  as  suggested  weakly  by  the  air  overpressure 
records  (Fig.  C.15 )•  Extension  of  integration  to  times  corresponding 
to  the  negative-phase  portion  of  the  stress  data  has  probably  had  lit¬ 
tle  effect  on  the  vertical  and  radial  velocity  and  displacement  curves, 
all  of  which  could  have  been  terminated  at  about  0.8  Bee  and  approxi¬ 
mately  zero  velocity  without  altering  the  early  portions  of  the  euxveB. 
Similar  correction  to  zero  velocity  either  at  the  positive  phase  or 
later  in  the  negative  phase  was  not  feasible  for  tangential  acceler¬ 
ometer  data,  partly  because  signal  strengths  were  20  per  cent  or  less 
of  Bet  range  after  about  0.6  sec. 

Arrival  of  the  main  shock  following  the  precursor  shock  at  about 
0.35  sec  is  evident  in  all  three  acceleration  components  for  Shot  10 
and  to  a  lesser  degree  in  the  derived  velocity  and  displacement  data. 
There  is  a  suggestion  of  a  similar  occurrence  in  the  Shot  1  data  shortly 
before  0.6  sec,  but  this  effect  is  not  so  well  corroborated  by  the 
stress  and  strain  curves  for  Shot  1  as  by  those  for  Shot  10, 

Peak  and  residual  displacements  seem  reasonable  with  the  possi¬ 
ble  exception  of  that  shown  fox'  the  upward  peak  of  EAV  displacement 
(0.82  in. )  for  Shot  1,  This  displacement,  greater  than  the  peak  down¬ 
ward  displacement,  may  have  resulted  from  incorrect  adjustment  of  the 
acceleration  curve.  The  adjustments,  which  were  of  the  order  of  0.03 
to  0.005  g,  were  very  small  but  often  affected  radically  later  parte 
of  the  displacement  curves. 


55 

SECRET- RESTRICTED  DATA 


The  most  optimistic  attempts  at  correlation  between  displace¬ 
ment-time  curves  and  the  earth  strain  curves  yielded  only  one  case  of 
apparently  similar  parameter-time  behavior.  The  vertical  displacement 
curve  (Fig.  C,2l)  and  the  vertical  strain  curves  (KS-l  and  ES-7)  for 
Shot  10  (Figs.  C.l6  and  C.l8)  are  of  generally  comparable  shape.  Per¬ 
haps  to  expect  correlation  of  these  two  parameters,  one  measured  by 
gages  which  are  very  sensitive  to  placement  conditions,  the  other  de¬ 
rived  by  a  double  integration  process  which  involved  arbitrary  zero 
adjustments,  is  ashing  too  much  of  the  data. 

Other  considerations  must  enter  into  speculation  on  the  rela¬ 
tion  between  strains  and  displacements,  however.  Strains  are  purely 
local  parameters  and  can  actually  represent  local  displacement  only  if 
the  variation  of  strains  over  a  volume  of  the  medium  1b  known.  Thus 
the  instantaneous  vertical  displacement  at  a  point  within  the  soil  rep¬ 
resents  the  integral  of  the  instantaneous  strains  induced  within  a 
column  of  soil  extending  from  the  point  of  measurement  to  some  depth 
below  which  loads  and  strains  become  vanishingly  small.  This  displace¬ 
ment  will  be  independent  of  the  strains  above  the  me&sui-ernent  point  at 
the  instant  under  consideration  and  will,  in  fact,  depend  on  conditions 
of  load  and  strain  in  the  overlying  material  only  as  these  conditions 
affect  and  have  affected  transmittal  of  load.  Correlation  of  displace¬ 
ment  with  strain  would  therefore  appear  to  require  knowledge  not  only 
of  tlie  strain  at  the  measuring  point  but  at  a  sufficient  number  of 
points  to  establish  strain  as  a  function  of  distance  and  the  Iocub  of 
vanishingly  small  strains. 

2.6  CONCLUSION 


Conclusions  are  divided  into  those  dealing  with  stress-strain 
characteristics  of  transiently  loaded  soil,  derivation  of  the  stress 
tensor,  and  derivation  of  velocity  and  displacement  information  from 
measured  accelerations. 

1.  It  is  feasible  to  install  existing  stress  and  strain  instru¬ 
ments  in  the  earth  at  depths  as  great  as  15  ft,  except  in  unusual  soil 
conditions . 

2.  Reproducibility  of  measurements  of  earth  stress  and  strain 
components,  using  the  methods  outlined  (as  expressed  by  the  average 
deviation  from  the  mean  of  duplicate  measurements),  is  better  than  16 
per  cent  for  stress  gages  and  about  25  per  cent  for  earth  strain  gages. 

3.  Stress-strain  data  from  component  measurements  in  the  earth 
provide  a  basis  for  estimating  energy  dissipation  caused  by  plastic 
defoliation  of  the  soil  under  transient  loading.  This  rate  of  energy 
loss  was  computed  to  be  of  the  order  of  30  pft-lb/cu  ft  at  the  5  ft 
depth  when  the  vertical  stress  is  of  the  order  of  20  to  30  psi  and 
300  pft-lb/cu  ft  when  the  vertical  stress  is  about  60  psi. 

4.  Earth  stress  tensors  and,  by  inference,  earth  strain  tensors 
in  the  form  of  principal  stresses  and  principal  strains,  may  be  derived 
from  measurements  of  properly  oriented  components  of  the  transient  earth 
stresses  and  earth  strains  produced  by  air  or  underground  explosions. 

The  assumption  of  instantaneous  equilibrium  of  stress  for  the  purpose 

of  using  analytical  procedures  based  on  static  stress  equilibrium  is 
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feasible  although  analysis  of  the  data  does  not  necessarily  indicate 
fulfillment  of  the  equilibrium  condition. 

5.  Use  of  earth  acceleration-time  data  for  derivation  of  velocity 
arid  displacement -time  information  is  feasible  within  certain  limits  of 
precision  of  the  data  and  duration  of  the  transient  phenomena,  as  es¬ 
tablished  by  several  previous  tests. M/  Correlation  between  earth 
displacement  data  so  derived  and  measured  earth  strains  is  not  clearly 
established  although  fair  agreement  was  indicated  in  one  case. 

It  is  interesting  to  consider  the  results  of  this  study  of  earth 
stress  and  strain  measurements  in  terms  of  the  applications  which  the 
techniques  might  find  in  solving  problems  of  underground  effects  of 
explosions  and  establishing  bases  for  theoretical  treatment  of  such 
problems.  Correlation  of  free -field  phenomena  with  damage  and  struc¬ 
tural  response,  is,  of  course,  necessary  to  usefulness  of  free-field 
measurements  and  is  not  yet  an  accomplished  fact,  possibly  because  of 
inadequate  knowledge  of  structural  response  and  of  the  coupling  of 
soil  and  structural  loads.  This  discussion  will,  however,  deal  with 
free-field  measurements  and  assume  their  usefulness  to  the  basic  prob¬ 
lems  of  damage.  The  assumption  will  also  be  made  that  methods  and 
techniques  of  installing  gages  in  the  earth  will  have  been  developed, 
at  least  to  such  a  point  that  realistic  corrections  can  be  made  in  the 
data  for  the  influence  of  perturbations  in  the  Btress  system  caused  by 
the  gages  and  by  the  disturbance  of  the  Boil. 

Air  bursts  of  nuclear  weapons,  particularly  those  at  heights  of 
the  order  of  one  or  two  fireball  radii,  induce  large  transient  stresses 
in  the  earth  within  a  stable  area  around  ground  zero.  Such  stresses 
are  evidently  predominantly  vertical  for  peak  incident  overpressures 
of  35  to  40  psi,  lateral  stresses  being  10  per  cent  or  less  of  the 
vertical  stresses.  There  is  no  reason  to  believe  that  this  situation 
would  change  at  higher  incident  overpressure  levels  except  possibly 
toward  an  even  lower  ratio  of  lateral  to  vertical,  stresses.  Conse¬ 
quently  free-field  measurements  concerned  with  damage  to  underground 
structures  from  air  burst  weapons  can  be  limited  to  the  vertical  com¬ 
ponent  of  various  phenomena  and  to  the  region  within  the  limits  of 
loading  which  may  produce  serious  structural  response.  The  principal 
problems  involve  the  rate  of  dissipation  of  the  energy  within  the  soil 
and  its  meaning  to  structural  response.  Correlation  of  vertical  btress - 
and  strain-time  data  as  a  function  of  depth  and  analysis  of  streBs-strain 
relations  as  a  function  of  depth  can  be  expected  to  give  useful  infor¬ 
mation.  An  additional  effect  suggested  by  Part  I  is  that  variation  of 
the  rate  of  stress  attenuation  or  energy  dissipation  with  depth  as  a 
function  of  peak  incident  overpressure  probably  merits  inclusion  in  a 
study  of  the  type  suggested  above.  The  apparent  increased  rate  of 
stress  attenuation  at  higher  incident  overpressures  might  well  be  a 
highly  significant  factor  in  the  effectiveness  of  earth  cover  protection 
for  underground  structures. 

All  these  effects,  to  be  useful  either  empirically  or  theoretically, 
must  be  correlated  with  structural  response  studies  of  the  type  initi¬ 
ated  by  Program  3.8  of  UPSHOT- KNOTHOLE. 5/ 

Application  of  free-field  measurements  of  earth  stress  and  strain 
to  underground  explosions  is  considerably  broader  than  for  air  bursts. 
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There  is  no  obvious  predominance  of  a  single  component  as  for  air  bursts. 
Radial  components  of  parameters  are  usually  greater  for  underground 
bursts,  but  vertical  components  are  frequently  nearly  as  -large  as  the 
radial  because  of  reflection  and  refraction  effects.  The  tangential 
component  may  be  negligible  in  seme  situations  so  that  radial  symmetry 
and  plane  stress  may  be  assumed.  However,  gross  inhomogeneities  of 
the  soil,  dissymmetry  and  development  of  shears  may  be  expected  to  make 
tangential  components  significant  in  many  cases.  Mo  sufficiently  com¬ 
prehensive  earth  stress -strain  measurements  have  been  made  in  connection 
with  underground  explosions  to  determine  satisfactorily  the  full  sig¬ 
nificance  of  the  various  component  measurements  to  an  understanding  of 
the  free -field  phenomena.  The  finite  size  of  gages  wad  the  number  re¬ 
quired  for  experimental  studies  of  stress  and  strain  tensor  measurement 
techniques  using  underground  bursts  of  moderately  small  underground 
high-explosives  charges  —  of  the  order  of  several  hundred  pounds  of 
TNT  —  require  gage  placement  at  such  relatively  remote  ground  ranges 
(in  terms  of  damage  criteria)  that  effects  of  soil  inhomogeneities, 
shears,  and  dissymmetry  may  be  obscured. 

Knowledge  of  the  stress  and  strain  tensors  as  a  function  of  posi¬ 
tion  with  respect  to  ground  zero  could  provide  a  more  realistic  basis 
for  estimates  of  underground  structural  damage  than  the  current  rule 
of  thumb.  Similarly  knowledge  of  the  stress-strain  pattern  would  be 
helpful  in  estimating  loading  on  structures.  These  are  free-field  meas¬ 
urements,  which,  to  have  usefulness,  must  be  correlated  with  structural 
response.  An  intermediate  step  in  this  correlation,  which  has  so  far 
been  neglected  and  which  is  probably  very  important,  is  measurement  of 
tensors  in  the  soil  immediately  adjacent  to  structures .  It  is  probably 
in  tills  zone  that  significant  differences  between  free-field  phenomena 
and  structural  loading  will  be  resolved.  The  plastic  nature  of  soil 
and  yielding  of  structural  elements  under  load  cause  a  serious  differ¬ 
ence  between  loading  and  response  of  soil  in  the  free-field  region  and 
the  corresponding  effects  at  a  structure.  Free-field  information  would 
also  serve  to  establish  boundary  values  for  theoretical  studies  of  dif¬ 
fraction  in  the  vicinity  of  underground  structures. 

It  does  not  seem  probable  that  a  sound  basis  for  design  or  damage 
estimates  for  underground  structures  will  be  available  until  the  stress- 
strain  reactions  of  ttie  soil  around  the  structure  are  understood  in  suf¬ 
ficient  detail  to  give  a  realistic  picture  of  the  loading  pattern  on 
the  structure. 


58 

SECRET  -  RESTRICTED  DATA 


APPENDIX  A 


AND  IMPULSE  INTENSITY-TIME  WAVE  FORKS  POP.  PROJECT  1.4a 


STRESS-TIME 


Fig.  A. 2  Earth  Stress  vs  Time  vs  Depth 
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Fig.  A 
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3  Earth  Stress  vs  Time  vs  Depth 
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Fig.  A. 5  Earth  Stress  vs  Time  vs  Depth 
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Fig.  A 


IA  |0  10  2  1  it  21  21 

r imr  |%r<  I 


6  Impulse  Intensity  vs  Time  vs  Depth 
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Fig.  A 
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9  Impulse  Intensity  vs  Time  vs  Depth 
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APPENDIX  B 


INSTALLATION  OF  GAGES  FOR  PROJECT  1.4b 


Earth  stress  and  earth  strain  gages  were  placed  in  trenches  5  ft 
deep  which  were  perpendicular  to  the  blast  lines  at  Stations  3-205  and 
F-28l.  Location  stakes  were  set  prior  to  excavation  to  show  the  center- 
line  of  the  trench  and  permit  reasonably  precise  positioning  of  the 
gages.  Each  trench,  as  excavated  by  a  backhoe,  had  a  bottom' width  of 
about  2  ft  and  a  width  of  about  '(  ft  at  the  surface .  The  length  of 
each  trench  was  sufficient  to  accommodate  the  full  complement  of  gages 
(Figs.  2.1  and  2.2) . 

The  bottom  of  each  trench  was  tamped  after  cutting  down  or  adding 
moist  backfill  as  needed  to  permit  installation  of  all  gages  with  re¬ 
sponse  axes  at  the  prescribed  depth  of  5  ft  within  +0.1  ft. 

The  position  of  each  gage  in  the  trenches  was  determined  from  two 
rows  of  location  stakes  (Fig.  B.l),  set  to  grade  at  3 -ft  intervals 
along  the  length  of  the  trenches,  one  row  7  ft  toward  ground  zero  and 
the  other  7  ft  away  from  ground  zero  with  respect  to  the  centerline. 

These  stakes,  undisturbed  by  excavation,  permitted  exact  positioning  of 
the  stress  gages  by  means  of  a  5  ft  plumb  line  dropped  from  the  mid¬ 
point  of  a  tape  stretched  between  corresponding  stakes  in  each  line . 
Strain  gages  were  placed  approximately  midway  between  adjacent  stress 
gages . 

Instruments  were  placed  in  the  trench  in  sequence  from  the  end 
nearest  the  blast  line.  Each  gage  was  oriented  with  respect  to  ground 
zero  by  means  of  a  sighting  tube  (Fig.  B.2)  and  relative  to  the  ver¬ 
tical  axis  by  means  of  a  spirit  level  and  protractor  (Fig.  B.3). 

Backfill  material,  consisting  of  soil  free  of  rock  and  gravel, 
was  mixed  with  water  and  allowed  to  season  for  about  a  day  to  approx¬ 
imately  optimum  compaction  consistency.  This  material  was  screened 
and  hand  tamped  around  each  instrument  (Fig.  B.H)  .  Backfill  was  placed 
in  thin  lifts,  about  ?.  in.  for  hand  tamping  and  6  to  12  In.  for  pneu¬ 
matic  tamping  (Fig.  B.5),  and  compacted  to  uniform  densities  somewhat 
greater  than  the  natural  density  of  the  surrounding  soil.  A  few  den¬ 
sity  determinations  of  the  undisturbed  soil  and  of  the  compacted  back¬ 
fill  were  made  for  control  purposes.  Backfill  was  placed  to  an  elevation 
of  about  2  ft  above  each  gage  before  the  next  one  was  installed. 
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Fig.  B.l  Instrument  Trench, 
Location  Stakes  and 
Gages  in  Prescribed 
Positions 


Fig.  B.2  Orientation  of  Earth 
Strain  Gage 


Fig.  B.3  Orientation  of  Earth  Fig.  B.4  Placement  of  Backfill 
Strain  Gage  Around  Gage 
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Earth  strain  gages  require 
more  meticulous  placement  than 
stress  gages.  The  very  short 
ranges  (0.010  in.  for  a  5  in. 
gage  length)  through  which  the 
gage  must  operate  after  comple¬ 
tion  of  all  backfill  and  the 
short  range  of  linear  response 
(0 .040  in.)  cf  the  transducers 
indicate  the  precision  with 
which  control  must  be  exercised 
during  placement.  Clamping 
devices  intended  to  restrain 
the  gages  in  proper  mechanical 
alignment  are  impractical. 

Requisite  control  of  strain 
gage  alignment  was  accomplished 
by  monitoring  the  output  of  each 
unit  during  its  placement  through 
a  portable  single -channel  car¬ 
rier-amplifier  unit  (Consoli¬ 
dated  Type  113B) -  The  amplifier 
output  was  under  observation  by 
the  man  compacting  the  backfill 
around  the  gage  throughout  the 
operation.  Alternation  of  tamping  operations  between  areas  which 
caused  extension  of  the  gage  and  those  which  produced  compression  per¬ 
mitted  thorough  compaction  and  attained  remarkably  good  ultimate  me¬ 
chanical  alignment.  This  monitoring  procedure  permitted  extension  or 
compression  of  u  gage  during  compaction  beyond  its  linear  range  with¬ 
out  loss  of  control. 

Pneumatic  tamping  replaced  hand  tamping  after  sufficient  fill  had 
been  placed  to  ensure  protection  of  the  gages  from  mechanical  damage. 
Soil  immediately  adjacent  to  the  gages  was  compacted  to  densities 
closely  approximating  that  of  the  surrounding  fill  and  exceeding  that 
of  the  undisturbed  soil. 

Considerably  greater  difficulty  was  experienced  in  control  of 
vertical  strain  gages  than  of  those  which  responded  to  other  compo¬ 
nents.  The  source  of  this  difficulty  appears  to  be  a  thin  layer  of 
soil  directly  beneath  the  upper  plate  of  the  vertically  placed  gage. 
Tills  soil  must  be  compacted  by  very  indirect  methods.  The  procedure 
found  most  effective  for  vertical  gages  was  to  compact  the  soil  be¬ 
tween  the  end  plates  sufficiently  to  extend  the  gage  through  several 
set  ranges  prior  to  final  compaction  of  the  overlying  soil.  This 
latter  operation  usually  followed  placement  of  a  little  less  than  1 
foot  of  fill  above  and  around  the  gage.  Compaction  operations  above 
the  gages  after  they  had  been  covered  by  more  than  1  foot  of  compacted 
fill  did  not  generally  result  in  permanent  mechanical  zero  changes  of 
more  than  a  few  tenths  of  set  range. 

The  two  vertical  earth  strain  gages  at  Sta  3-205  were  not  very 
well  zeroed.  The  final  set  of  ES-1  was  off  true  mechanical  zero  by  an 
appreciable  fraction  of  a  set  range,  and  that  of  ES-Y  was  even  more 


Fig.  B.5  Tamping  Backfill 
Over  a  Gage 
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seriously  out  of  alignment.  However,  neither  departure  exceeded  the 
linear  range  of  the  gage.  The  corresponding  gages  at  Sta  F-28l  were 
zeroed  so  that  the  final  set  was  within  a  very  small  fraction  of  a  set 
range  from  true  mechanical  zero. 

Compaction  of  the  final  3  ft  of  fill  above  the  gages  did  not  af¬ 
fect  the  mechanical  setting  of  the  strain  gages.  All  gages  responded 
to  the  tamper  loading  but  returned  to  their  initial  as-placed  zeros. 
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APPENDIX  C 


PARAMETER  TIME  RECORDS 


Fig.  C.l  Earth  Stress,  Strain,  and  Air  Overpressure  vs  Time  -  Shot  1 
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Fig,  C.2  Earth  Stress  and  Strain  vs  Time  -  Shot  1 
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Fig.  C.3  Earth  Stress  and  Strain  vs  Time  -  Shot  1 
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Fig.  C 


,4  Earth  Stress  and  Strain  vs  Time  - 
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Shot  1 
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Fig.  C.5  Vertical  Ground  Acceleration,  Velocity,  and  Displacement 
vs  Time  -  Shot  1 
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Fig.  C 
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7  Tangential  Ground  Acceleration,  Velocity,  and  Displacement 
vs  Time  -  Shot  1 
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g.  C.8  Air  Overpressure  and  Earth  Stress  vs  Time  —  Shot  9 
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Fig.  C.9  Earth  Stress  and  Strain  vs  Time  -  Shot  9 
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Fig.  G.10  Earth  Stress  end  Strain  vs  Time;  -  Shot  9 
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Vertical  Ground  Acceleration,  Velocity,  and  Displacement 
Time  -  Shot  9 
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C.13  Radial  Ground  Acceleration,  Velocity,  and  Displacement  vs 
Time  -  Shot  9 
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Fig.  C.14  Tangential  Ground  Acceleration,  Velocity,  and  Displacement 
vs  Time  -  Shot  9 
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Fig.  C.15  Air  Overpressure  and  Earth  Stress  vs  Time  -  Shot  10 
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Fig.  C. 
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]6  Earth  Stress  and- Strain  vs  Time  -  Shot  10 
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Fig.  C.17  Earth  Stress  and  Strain  vs  Time  -  Shot  10 
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Fig.  C.18  Earth  Stress  and  Strain  vs  Time  -  Shot  10 
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C.21  Vertical  Ground  Acceleration,  Velocity,  aftd  Displacement  vs 
Time  -  Shot  10 
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Fig.  C.22  Radial  Ground  Acceleration,  Velocity,  and  Displacement  vs 
Time  -  Shot  10 
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APPENDIX  D 


STRESS  TENSOR  -  PRINCIPAL  STRESS  DERIVATION 


Stress  is  the  normal  or  tangential  component  of  pressure  acting 
on  a  surface  within  a  solid  or  quasi -solid.  Stresses,  like  pressures, 
represent  intensity  of  force  per  unit  area.  The  stress  tensor  is  the 
mathematical  expression  which  represents  completely  the  condition  of 
stress  acting  through  a  point. 

The  stress  tensor  is  defined  by  all  stresses  acting  on  three 
mutually  orthogonal  planes.  In  the  general  case  this  definition  re¬ 
quires  knowledge  of  nine  stress  components,  one  normal  to  each  plane 
and  two,  usually  mutually  perpendicular,  tangent  to  each  plane.  How¬ 
ever,  equilibrium  of  stresses  is  necessary  for  analysis  of  the  system, 
and  this  condition  requires  equality  of  three  pairs  of  tangential 
stresses.  Consequently  the  general  stress  tensor  is  defined  by  six 
stress  components.  There  is,  in  addition  to  this  definition  of  the 
tensor,  a  special  case  which  results  from  the  fact  that  for  any  point 
in  a  solid  in  which  stresses  are  in  equilibrium  there  exists  one  orien¬ 
tation  of  three  mutually  orthogonal  planes  such  that  all  tangential 
stresses  on  them  vanish.  These  are  the  principal  planes.  This  special 
case  completely  defines  the  stress  tensor  in  terms  of  only  the  three 
normal  stresses  on  the  principal  planes,  i.e.,  the  principal  stresses. 

The  purpose  of  this  appendix  is  to  show  how  data  from  practical 
field  measurements  of  earth  stress  components  may  be  treated  mathemati¬ 
cally  to  define  the  stress  tensor  in  terras  of  the  three  principal 
stresses.  Figure  D.l  defines  the  terminology  and  geometry  of  the 
stresses  involved.  As  noted  in  Part  II,  the  direction  of  a  plune  is 
defined  by  the  normal  to  it;  thu3  direction  of  a  plane  and  its  normal 
stress  ere  defined  by  the  some  symbol.  Three:  mutually  orthogonal 
planes  through  the  point,  0,  of  the  diagram,  are:  the  V-plane  which 
is  horizontal  and  below  the  ground  surface  at  the  depth  of  the  instru¬ 
ments  as  depicted  in  the  upper  part  of  the  diagram;  the  R-plane  which 
is  vertical  and  normal  to  the  radius  from  ground  zero  at  0;  and  the 
T-plane  which  is  also  vertical  and  normal  to  the  tangential  axis 
through  0.  These  planes  are  chosen  for  discussion  since  the  normal 
stresses  on  them,  by,  dj{,  and  dip,  at  0  correspond  to  three  of  the 
measured  stresses  discussed  in  Section  2.5.3*  The  upper  portion  of 
the  diagram  shows  a  fourth  plane,  which  for  purposes  of  analysis 
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STATION  F-  Jil 


represents  a  principal  plane,  although,  its  orientation,  as  shown,  is 
wholly  arbitrary  and  was  chosen  only  to  ensure  contrast  with  the  V-, 

R-,  and  T-planes. 

The  tetrahedron  in  the  lower  part  of  Fig.  D-l  represents  an  in¬ 
finitesimal  portion  of  the  four  planes  at  0.  The  three  stresses  on 
each  of  the  orthogonal  planes  at  0  are  represented  by  arrows,  and  the 
normal  stress  on  the  arbitrary  plane  is  represented  by  a  perspective 
arrow.  The  stress  arrows  in  the  V-,  R- ,  and  T-planes  were  drawn  as 
tension  for  clarity  in  the  drawing  but  should  all  point  inward  to  rep¬ 
resent  compression  as  does  the  dp-arrow.  Compressive  stresses  are 
taken  positive  in  accordance  with  soil  mechanics  convention  and  con¬ 
trary  to  general  stress  analysis  usage.  Normal  stress  is  identified 
by  o  with  a  subscript  representing  the  plane  on  which  they  act;  shear 
stress  is  designated  by  a  with  two  subscripts,  indicating  the  direc¬ 
tion  of  the  shear  in  terms  of  (l)  the  axis  to  which  it  is  parallel  and 
(2)  the  plane  in  which  it  acts.  Thus  dy  is  the  normal  stress  on  the 
V-plane  and  -tpy  is  the  shear  parallel  to  the  R-axis  in  the  V-plane. 
op  is  the  normal  stress  on  the  arbitrary  or  P-plane,  which,  because  of 
its  definition  as  a  principal  plane,  has  no  shear  stresses  associated 
with  it. 

The  stress  tensor  is  defined  by  the  nine  components  ay,  dp,  d<p> 
tpy,  -vpy j  -typ*  iTR>  and-  ^RT*  However,  since  stress  equilibrium 

must  be  assumed  for  analysis,  it  is  necessary  that  ayp  =  -tRy ,  typ  =  Trjy , 
and  tprji  =  trj^;  whence  the  stress  tensor  is  defined  by  six  components, 
dy,  Or,  d>j,  iyp,  Typ»  and 

Equilibrium  of  stresses  at  0  requires  that  all  forces  on  the  in¬ 
finitesimal  tetrahedron  be  in  equilibrium.  These  forces  are  the  pro¬ 
duct  of  the  stresses  by  the  areas  over  which  they  act.  If  A  represents 
the  triangular  area  of  the  P-plane  included  as  a  face  of  the  tetrahedron, 
then  Ak,  Am,  and  An  are  the  triangular  areas  of  the  faces  In  the  V-,  R-, 
and  T-planes  and  k,  m,  and  n  are  directional  cosines  of  the  P-plane  with 
respect  to  the  V-,  R-,  and  T-axes,  i.e.,  k  =  cob  (P,V),  m  -  cos  (P,  R), 
and  n  =  cos  ^P,T).  The  force  Adp  on  the  P-plane  may  be  resolved  into 
components  parallel  to  the  axes,  and  each  of  these  component  forces  must 
be  equal  in  magnitude  and  opposite  in  direction  to  the  corresponding 
forces  derived  from  the  stresses  on  the  V-,  R-,  and  T-faces  of  the  tetra¬ 
hedron.  Thus  for  equilibrium  three  equations  of  the  type 

dp  Ak  -  dy  Ak  -  -typ  Am  -  An  =  0 

must  be  satisfied.  The  equilibrium  equations  from  which  the  area  terms 
have  been  cancelled  :ire 


(Op  -  <V 

)  k  - 

W  -  V  =  ° 

(D.l) 

-^VR  k  + 

(<»p 

-  «R>  m  "  W  =  ° 

(D.2) 

-‘‘VT  k  ' 

XRT 

m  +  (dp  -  dT)  n  =  0. 

(D.3) 
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This  group  of  simultaneous  equations  may  he  solved  to  yield  the 
directional  cosines  if  the  values  of  tip  are  known.  These  solutions 
can  differ  from  zero  only  if  their  determinant  vanishes.  That  is, 
only  if 

(°P  -  ^vr  "Vr 

“tVR  (°P  "  °r)  -tTR  =  °* 

~tVT  -,cRT  °T^ 

This  determinant  may  be  restated  as  the  cubic  equation  in  Op, 

3  2 

dp  +  Pdp  +  qdp  +  r  =  0, 

where 

p  =  -  (dy  +  dR  +  dT) 

2  2  2 

a  =  +  (“Vr  +  dR°T  +  dVdT  "  'CVR  '  ’'rT  “  TVr^ 

r  =  -  (dydpd^,  +  2TVRXRTtVT  "  °VCRT  "  °RTVT  '  VW' 

The  three  roots  of  this  equation,  d^,  dg,  and  d^,  are  normal  stresses 
on  three  mutually  orthogonal  planes  in  which  Bhear  stresses  vanish  and 
arc  consequently  the  principal  stresses  at  0. 

Solution  of  the  cubic  equation  requires  elimination  of  the  square 
term  and  conversion  of  the  resultant  to  a  quadratic  equation  in  cubes. 
Thus  if  dp  =  x  +  c,  Eq  D.5  becomes 

*3  2  2  o  p 

x”5  +  (3c  +  p)  x  +  (3c  +  2pc  +  q)  x  +  (cJ  +  pc  +  qc  +  r)  =  0, 
and  if  3e  +  p  =  O,  the  x  -term  vanishes  and  c  a  -p/3*  Eq  D-5  is  now 

x^  +  ax  +  b  =  0.  (0.6) 

(D.7) 

(D.8) 
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(D.4) 


(D.5) 


Equation  D.6  may  tie  converted  to  e.  quadratic  form  by  placing  x  =  y  +  z, 
whence 

y^  +  +  (3yz  +  a)  (y  +  z)  +  b  =  0  (D.9) 


from  which  the  (y  +  z)-term  may  be  eliminated  by  making  (3yt  +  a)  vanish. 
That  is,  if  3yz  +  a  -  0,  z  -  -  a/3y,  Sq  D.9  becomes 


y 


3 


b  =  0 


or 


6 

y 


+  by~ 


(D.10) 


The  last  equation  is  a  quadratic  in  y  which  has  the  more  familiar 
form  s2  +  bs  -  a3/27  =  0  if  y3  is  replaced  by  s . 

The  solution  of  this  quadratic  equation  is 


and  similarly, 


(D.ll) 


(D.12) 


The  solution  of  Eq  D.6  is  the  sum  of  the  three  cube  roots  of  Eqs 
D.ll  and  D.12  since  x  =  y  +  z.  However,  there  will  be  three  possible 
sets  of  roots,  depending  on  whether  the  quantity  beneath  the  radical 
in  !3qs  D.ll  and  D.12  is  greater  than,  equal  to,  or  less  than  zero.  For 
if  b2/4  +  a3/27  >  0,  only  one  root  will  be  real  and  the  other  two  con¬ 
jugate  complex;  if  b^/4  +  a.3/27  --  0,  nil  roots  will  be  real,  but  two 
will  be  equal;  and  if  b2/4  +  n3/27  <  0,  all  roots  will  be  real  and  un¬ 
equal.  The  last  of  these  conditions  is  the  only  one  that  is  of  real 
concern  to  the  derivation  of  principal  stresses.  Unfortunately  this 
case  is  also  the  one  which  is  algebraically  irreducible;  algebraic 
manipulation  of  the  equation  lends  only  to  identical  equations  because 
of  the  conjugate  complex  character  of  Eqs  D.ll  and  D.12  when  b2/4  + 
3.3/27  is  less  than  zero.  It  is  feasible,  however,  to  derive  trigono¬ 
metric  solutions  for  Eqs  D.ll  and  D-12  which  yield  usable  information 
for  evaluation  of  x  and  subsequently  of  Op.  2 

The  quantities  y3  and  z3  are  conjugate  complex  functions  if  b  /4 
+  a3/27  <  0,  i.e.,  if  ai/27  is  negative  and  numerically  greater  than 
b^/4  when  the  latter  is  positive.  Equations  D.ll  and  D.12  may  be 
rewritten  in  the  trigonometric  form 


=  w  (cos  9  +  i  sin  ©)  (D.13) 
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,3 

where,  If 


(cos  0-1  sin  O) 


u 


w  cos  a  = 


-b 

2 


and 

v  =  w  Bin  6  -  V"  V*”  +  ff)  '  (D*i5) 

and 

cos  a  -=  ^  .  (d.i6) 

Equations  D.13  and  D.lU  may  now  be  evaluated  by  taking  the  cube 
roots  of  the  trigonometric  expressions 

■  w1^  (coa  9  +  i  sin  9)X/3  x 

=  w1/3  (cos  9  -  i  Bin  O)1^  x  (i)^/3  , 

where  (l)1^  has  the  values  1,  (cos  2n/3  +  i  sin  2w/3)>  and  (cos  4w/3 

+  i  sin  Un/3).  The  result  is  three  values  each  for  y  and  z: 

yl  = 
y2  = 


Z1  = 


Z3 

100 


(=0,  | 

uV3  £ 


cos 


Icos 


+  1  sin 

9  +  2r 


f) 


3 

9  +  4ir 


+  i  sin 
+  i  sin 


v1^  L 

Lb  §  -  1  Bln  f) 

1/3  f  9  +  2n  ,  , 

w  '  J  Icos  — - — u  -  i  sin 


9  4-  2ir 
3 

e  +  Uw1 


3 

9  •*■  4w 
3 


-  i  sin 


+  2u\ 

3  J 

+  4n\ 

3  J 
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Consequently  the  solutions  of  Eq  D.6  are  the  sums  of  these  ex¬ 
pressions  : 


X^  —  c.  w 


..1/3 


0 

cos  •?  , 


„2  -  2  v1^  oos  ii-Si  -  -2  »1/3  cos 


*3  -  2  v1*  cos  S-Sil  «  ,2  ul/3 


COS 


(«*  •  f)  ■ 

(60°  +  !)  • 


*2  1  /p 

However,  since  w  =  {-er /2lY'  and  dp=x+c=x-  p/3,  the 
solutions  for  Eq  D.‘>  are 


’i  ■ 2  (f)1/2  “s  (!)  - 


»2*-2 


(f)1/2  -  ka°  ■  !)• 


(D.17) 

(D.18) 


and 


(f)L/2  =os  (60°  +  "  3  *  (D.19) 


Equations  D.17,  D.l8,  and  D.19  are  the  three  principal  stresses  at 
the  origin  0  of  Fig.  D.l.  They  may  be  evaluated  completely  from  meas¬ 
ured  values  of  the  normal  stre36,  dy,  Op,  and  dT,  and  measured  or  de¬ 
rived  values  of  the  shear  stresses  iyp,  Tpip,  and  ty^  at  that  point 
through  the  definition  of  the  quantities  p,  q,  and  r  for  Eq  0.9,  of  a 
and  b  in  Eqs  D.7  and  D.8,  and  the  expression  for  cos  G  in  Eq  D.lb.  The 
subscripts  used  to  identify  the  three  stresses  of  Eqs  D.17,  D.18,  and 
D.19  are  those  customarily  used  to  distinguish  the  major,  intermediate, 
and  minor  principal  stresses,  a  distinction  based  upon  magnitude.  The 
relative  magnitudes  will  depend  upon  both  the  magnitudes  and  the  signs 
of  the  observed  normal  and  shear  stresses  and  will  not  necessarily  re¬ 
sult  in  values  such  that  dp  >  dg  >  dg  in  the  notation  given  above. 
Consequently  it  may  be  convenient  in  practice  to  interchange  the  sub¬ 
scripts  of  the  principal  stresses  given  by  Eqs  D.17,  D.l8,  and  D.19  to 
conform  with  customary  usage. 

A  further  problem  in  application  of  the  analytical  methods  to  data 
from  earth  stress  measurements  is  that  introduced  by  the  lack  of  suit¬ 
able  instruments  for  direct  observation  of  transient  earth  shear 
stresses  (Section  2.5.3).  This  problem  has  been  solved  by  utilizing 
the  normal  stresses  measured  on  two  orthogonal  planes  to  derive  the 
shear  stresses  on  the  plane  which  bisects  the  angle  between  them.  Thus 
normal  stresses  dyp  and  dy^  DA,  measured  on  planes  directed  at  45° 
downward  toward  ground  zero  ond  downward  away  from  ground  zero  between 
the  vertical  and  radial  planes,  were  used  to  derive  the  shear  stresses, 
Typ  and  -tRy,  on  the  radial  and  vertical  planes  by  means  of  the  equation 
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'VR  “  ‘‘RV 


(D.20) 


vr  nr 


°VR  DA 

2 


The  final  step  in  defining  the  stress  tensor  by  the  principal 
stresses  is  determination  of  the  direction  of  these  stresses  with  re¬ 
spect  to  the  coordinate  axes .  These  directions  are  expressed  as  three 
directional  cosines  for  each  stress,  representing  the  angles  between 
the  stress  and  each  axis.  A  directional  cosine  is  the  ratio  of  a  di¬ 
rectional  number  to  the  unit  vector  representing  the  principal  stress, 
and  directional  number  is  the  projection  of  the  unit  principal  stress 
vector  on  a  coordinate  axis. 

The  factors  k,  m,  and  n  of  the  simultaneous  Eqs  D.l,  D.2,  and  D-3 
were  defined  as  the  directional  cosines  of  the  P-plane.  Solutions  of 
these  equations  will  be  directional  numbers _whioh  by  normalization  are 
converted  to  directional  cosines.  Thus  if  k,  m,  and  n  are  solutions  of 
the  simultaneous  equations,  the  directional  cosines  are  given  by 


.  _  _ k 

K  ~  7,2  ^  -2"  '  -2V./2 

Ik  +  m  +  n  j  ' 

_ m 

m  _  ....  _2-  ^  _2v1y2 

Ik  +  m  +  n  V 

and  '  ’ 

n  =  -rg  Ig  _2Vl/2  • 

|k  +  m  +  n  IT' 

The  determinant  A  (Eq  D.4)  of  the  simultaneous  equations  must 
vanish  if  their  solutions  are  different  from  zero.  The  minors  of  the 
three  terms  of  any  row  of  the  determinant  can  be  shown  to  be  solutions 
of  the  simultaneous  equations  if  they  satisfy  the  equations.  Thus  if 
the  minors  of  the  terms  of  the  first  row  (op  -  by)*  ".^VR*  "  ^VT 

ore  k,  in,  and  n,  then  from  Eq  D.l  substitution  of  k,  in,  and  n  for  k,  m, 
and  n  gives  (tip  -  <jy)  k  -  -typin  -  typin  which  must  vanish  since  it  is 
exactly  the  determinant  A .  Similar  substitutions  in  Eqs  D.2  and  D.3 
give 


both  of  which  are  determinants  in  which  two  rows  contain  identical 
terms  and  must  consequently  vanish.  Thus  the  minors  of  the  first  row 


(D.21) 

(D.22) 

(0-23) 
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satisfy  the  simultaneous  equations  and _ are  solutions.  The  directional 
numbers  for  the  P-plane  are  therefore  it,  5,  and  n,  the  minors  of  the 
first-row  terms  of  the  determinant,  and  corresponding  directional  co- 
<  sines  are  given  by  Eqs  D.21,  D.22,  and  D.23. 

The  P-plane  has  been  defined  as  one  of  the  principal  planes  and 
Cp  as  the  corresponding  principal  stress.  Therefore,  the  directional 
r  cosines  of  the  principal  stress  defining  the  angles  between  it  and  the 

coordinate  axes  V,  H,  and  T  are 

k  = 


m  = 

and 


n  = 

where  the  subscript  P  has  the  values  ],  ?,  and  3  as  the  major,  inter¬ 
mediate,  and  minor  principal  stresses  are  used  in  determining  the  di¬ 
rectional  numbers . 
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